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Abstract 

Evidence of aerosol intercontinental transport (ICT) is both widespread and compelling. 
Model simulations suggest that ICT could significantly affect regional air quality and 
climate, but the broad inter-model spread of results underscores a need of constraining 
model simulations with measurements. Satellites have inherent advantages over in situ 
measurements to characterize aerosol ICT, because of their spatial and temporal 
coverage. Significant progress in satellite remote sensing of aerosol properties during the 
Earth Observing System (EOS) era offers opportunity to increase quantitative 
characterization and estimates of aerosol ICT, beyond the capability of pre-EOS era 
satellites that could only qualitatively track aerosol plumes. EOS satellites also observe 
emission strengths and injection heights of some aerosols, aerosol precursors, and 
aerosol-related gases, which can help characterize aerosol ICT. After an overview of 
these advances, we review how the current generation of satellite measurements have 
been used to (1) characterize the evolution of aerosol plumes (e.g., both horizontal and 
vertical transport, and properties) on an episodic basis, (2) understand the seasonal and 
inter-annual variations of aerosol ICT and their control factors, (3) estimate the export 
and import fluxes of aerosols, and (4) evaluate and constrain model simulations. 
Substantial effort is needed to further explore an integrated approach using measurements 
from on-orbit satellites (e.g., A-Train synergy) for observational characterization and 
model constraint of aerosol intercontinental transport and to develop advanced sensors 
for future missions. 
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1. Introduction 

Atmospheric aerosols are solid or liquid particles suspended in the atmosphere, exhibiting 
great variability not only in concentration but also in size, composition, and shape. They 
are emitted directly into the atmosphere as primary aerosols, or are secondary aerosols, 
formed through chemical transformations from a variety of natural and man-made 
gaseous emissions. In high concentrations, aerosol, also known as particulate matter or 
PM, can be hannful to human health (Pope et al., 2002). Aerosols affect weather and 
climate by scattering and absorbing solar radiation (so-called “aerosol direct effects”) 
(McCormick and Ludwig, 1967) and by modifying cloud properties, amount, and 
evolution (collectively referred to as “aerosol indirect effects”) (Twomey, 1977; Gunn 
and Phillips, 1957; Albrecht, 1989). Absorption of solar radiation by particles changes 
the atmospheric stability and reduces the surface fluxes, which can change atmospheric 
circulations (Ramanathan et al., 2005; Lau and Kim, 2006; Lau et al., 2009; Zhang et al., 
2009). Dust aerosol rich in iron also affects climate indirectly by changing 
biogeochemistry (Jickells et al., 2005). Aerosol also impairs visibility (White, 1976) and 
affects crop yields (Chameides et al., 1999). 

Aerosols are removed from the atmosphere through gravitational settling, dry deposition, 
and wet scavenging, processes that typically produce lifetimes of hours to days in the 
atmospheric boundary layer and of weeks in the upper troposphere. Thus aerosols emitted 
from one continent are often transported long distances to another, in particular when 
aerosols are pumped out of the boundary layer. There is mounting evidence for 
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intercontinental and even hemispheric transport, provided by long-term surface 
monitoring networks, in-situ measurements from intensive field campaigns, and 
especially satellite observation, backed by model simulations. It was determined that the 
recurring phenomenon of “Arctic haze” in late winter and early spring is associated 
mainly with long-range transport of pollution and dust sources in Europe and Asia (Rahn 
et ah, 1977; Shaw et ah, 1983). Saharan dust routinely reaches the Amazon, the 
Caribbean Sea and southeastern U.S. after sweeping over the Atlantic Ocean (; Swap et 
ah, 1992; Prospero, 1999; Prospero et ah, 2005). Asian dust and pollution have long been 
observed in the western U.S. and beyond (Jaffe et ah, 1999; Biscaye et ah, 2000; 
VanCuren and Cahill, 2002; Grousset et ah, 2003; Stohl et ah, 2007). Because of the 
long-range transport, aerosols emitted or formed in one region could have important 
implications for climate, weather, and air quality in downwind regions/continents. As 
such, the long-range transport of air pollution has been an active topic of scientific 
research for several decades and has also been debated for formulating environmental 
policies and treaties (Holloway et ah, 2003). One of the most recent coordinated efforts is 
the establishment of the Task Force on Hemispheric Transport of Air Pollutants (HTAP) 
by the United Nations Economic Commission for Europe, aiming to understand the 
growing body of scientific evidence of intercontinental transport and assess its impacts on 
air quality, climate, and ecosystems (HTAP, 2010). 

Impacts of aerosol intercontinental transport (ICT) on regional air quality and climate 
change are potentially large, as estimated by chemical transport models (CTM) (HTAP, 
2010); however, large uncertainties are associated with the ensemble CTMs employed in 
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the assessment. Figure 1 shows the relative annual intercontinental response (RAIR) for 
surface aerosol concentration, a quantity related to air quality, and aerosol optical depth 
(AOD), a quantity related to aerosol climate forcing. RAIR is calculated from the source- 
receptor relationship experiments of multiple global aerosol models (HTAP, 2010; Yu et 
al., 2012a). By definition, RAIR in a receptor region represents the percentage 
contribution of the intercontinental transport of foreign emissions relative to the sum of 
foreign and domestic emissions (HTAP, 2010). Clearly, the aerosol import from ICT 
makes significant contribution to climate forcing and air quality. RAIR for surface 
concentration ranges from 5-20%, whereas RAIR for AOD is larger, having a range of 
10-30%. Values of RAIR depend on both region and component. On the other hand, 
model simulations show a large spread, as reflected by the calculated standard deviations 
of RAIR and represented by the error bars in Figure 1. The large model diversity is likely 
associated with differences in the model treatment of numerous aerosol processes, as 
extensively documented in a number of inter-comparisons between models and 
observations (e.g., Kinne et al., 2006; Schulz et al., 2006; Textor et al., 2006; Textor et 
al., 2007; Shindell et al., 2008; Koch et al., 2009; Prospero et al., 2010). 

There is a pressing need for reliable observations that constrain and improve model 
assessments of aerosol ICT and its regional impacts. Such observations need to meet 
some primary requirements. Measurements should have adequate spatial and temporal 
coverage to allow for tracking the evolution of aerosol transport over a global or 
hemispherical scale with a daily or hourly frequency. Measurements also should provide 
or can be used to derive highly accurate, quantitative information about aerosol 
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horizontal and vertical distributions, such as mass concentration, mass loading, aerosol 
optical depth (AOD), and vertical profiles of aerosol extinction. In the case of remote 
sensing observations, measurements should distinguish aerosol types, which is essential 
for the conversion between aerosol optical measurement (e.g., AOD and extinction or 
backscatter coefficient) and aerosol mass loading or concentration. Aerosol mass 
extinction efficiency and its dependence on relative humidity depend on aerosol type. 
Different types of aerosol also tend to follow different transport paths and exert different 
impacts on climate and human health. 

Large-scale satellite measurements, complemented by regional and local aerosol 
characterizations from surface stations and field campaigns, can meet these fundamental 
requirements. Surface stations can offer detailed albeit local measurements of surface 
concentrations, AOD (Holben et ah, 1998), and chemical (Malm et al., 2004), 
microphysical, and radiative (Dubovik et al., 2002) aerosol properties. Intensive field 
experiments can provide comprehensive snapshots of regional aerosols through 
implementing coordinated and multiplatform observation strategies (Russell et al., 1999; 
Huebert et al., 2003; Jacob et al., 2003; Reid et al., 2003; Tanre et al., 2003; Li et al., 
2007; Singh et al., 2009; Jacob et al., 2010). For all of their advantages, field campaigns 
and surface stations are inherently limited by the relatively short duration and/or small 
spatial coverage, which alone are not adequate for assessing the temporal and spatial 
variations of cross-ocean aerosol transport. Even long-term surface networks, such as the 
Aerosol Robotic Network (AERONET) (Holben et al., 1998), are substantially limited in 
geographical coverage, particularly over oceans where intercontinental transport occurs. 
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This lack of oceanic coverage is being addressed in part by the Marine Aerosol Network 
(MAN) (Smirnov et ah, 2009); but even so, marine aerosol coverage, both in space and 
time, remains exceedingly sparse. Satellite aerosol remote sensing can augment field 
campaigns and surface networks by expanding temporal and spatial scales because of the 
inherent advantage of daily or hourly measurements with global or continental coverage 
over decadal-scale durations. Indeed, ever since meteorological satellites were put into 
orbit, satellite imagers have been used to detect the aerial extent and motion of large- 
scale aerosol plumes, primarily in a qualitative or at most semi-quantitative way (e.g., 
Fraser, 1976; Lyons et al., 1978; Fraser et al., 1984; Chung et al., 1986; Ferrare et al., 
1990; Dulac et al., 1992). In recent decades and particularly during the Earth Observing 
System (EOS) era, substantial progresses have been achieved in satellite remote sensing 
of aerosols (see summaries by Yu et al., 2009a; Kahn, 2011; and references therein). 
Satellites are currently providing global land and ocean measurements of AOD with 
much improved quality. Global measurements of aerosol vertical distributions in cloud- 
free atmosphere and above cloud are also emerging (Winker et al., 2010). Enhanced 
remote sensing capabilities for constraining aerosol shape, size, and absorption allow for 
the characterization of aerosol type (Tanre et al., 2001; Higurashi and Nakajima, 2002; 
Kaufman et al., 2005a; Kahn et al., 2007; Tanre et al., 2011). Several satellite sensors are 
also providing value-added observations for aerosol ICT by measuring trace gases such 
as SO 2 and CO, that are aerosol precursors or originate from the same sources as 
particles, and providing means for identifying major aerosol source regions and 
estimating emissions (e.g., Veefkind et al., 2011; Ichoku et al., 2012). These advances 
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have provided an unprecedented opportunity of going from qualitative tracking to 
increasingly quantitative characterization and estimates of aerosol ICT. 

In this paper, we review: (a) the advances in satellite aerosol measurement over the last 
decade (section 2), (b) how the advanced satellite measurements have been used to 
improve the understanding of aerosol ICT variations on daily, seasonal, and decadal 
timescales and to estimate aerosol exports and imports (section 3. 1-3.3), and (c) the 
progress of using satellite measurements to evaluate and constrain model simulations 
(section 3.4). Major conclusions and recommendations are summarized in section 4. 

2. Recent advances in satellite aerosol measurements 

Satellite imagers have been providing a wealth of qualitative or semi-quantitative 
information for intercontinental and hemispheric aerosol transport since the 1970s (e.g., 
Prospero and Carlson, 1972; Lyons et al., 1978; Chung, 1986; Herman J. et al., 1997). 
Some studies also attempted to estimate columnar atmospheric density and transport of 
desert dust (Fraser, 1976; Mekler et al., 1977; Sirocko and Sarnthein, 1989; Dulac et al., 
1992), sulfate pollution (Fraser et al., 1984), and fire smoke (Ferrare et al., 1990) using 
early satellite measurements. However, such estimates were subject to large uncertainties 
because of poor data accuracy and the lack of constraints on aerosol microphysical 
properties and vertical distributions. Thus, quantitative assessments of aerosol 
intercontinental transport became possible only recently as a result of the much improved 
measurement accuracy and enhanced measurement capabilities achieved mainly in the 
EOS-era. Table 1 lists major satellite observables currently available for studying aerosol 
ICT. In the following we provide an overview of major advances achieved primarily 
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over the last decade, including improved accuracy and coverage for AOD, the emerging 
climatology of aerosol vertical distribution, enhanced capabilities for characterizing 
aerosol properties, and the feasibility of synergistic aerosol characterization. 

2.1. Improved accuracy and coverage for AOD measurements 

AOD is the quantity most often provided by satellite remote sensing. It is a measure of 
the integrated aerosol extinction (scattering and absorption) through the atmosphere, 
which is related to aerosol mass loading, dependent upon aerosol composition and 
meteorological conditions. Early aerosol monitoring from space used data from sensors 
that were designed for other purposes, e.g., the Advanced Very High Resolution 
Radiometer (AVHRR), the Total Ozone Mapping Spectrometer and (TOMS), and the 
Geostationary Operational Environmental Satellite (GOES). They have provided a multi- 
decadal climatology of AOD that has significantly advanced the understanding of aerosol 
distributions (e.g., Herman J. et al., 1997; Husar et al., 1997; Mishchenko et al., 1999; 
2003; Geogdzhayev et al., 2002; Torres et al., 2002; Zhao et al., 2008). However AOD 
accuracy from these early sensors is generally low. Major sources of AOD uncertainty 
are associated with instrument calibration, band characteristics, spatial resolution, cloud 
contamination, surface reflectance, and assumed aerosol optical models (Zhao et al., 
2008). 

In the past decade, significant progress has been made in reducing these uncertainties. 
On-orbit calibration has been implemented widely. Improvements in pixel resolution and 
spectral range make the cloud screening more robust than before (Martins et al., 2002), 
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though cloud contamination remains a major issue for aerosol remote sensing (Kaufman 
et al., 2005b; Zhang et al., 2005; Kahn et al., 2009). Several techniques based on 
enhanced sensor capabilities such as new channels (Kaufman et al., 1997; Hsu et al., 
2004), multi-angle views (Diner et al., 1998; Martonchik et al., 1998; 2002; 2009) and 
polarization (Hennan M. et al., 1997; Deuze et al., 2001), have been developed to 
improve the separation of land surface from atmospheric reflectance and hence, the 
accuracy and coverage of AOD retrievals. The advanced datasets of land surface 
properties have also been used to better characterize and quantify surface contributions in 
the aerosol retrieval algorithms (Li et al., 2005; Lyapustin et al., 2011). The accumulation 
of detailed aerosol property measurements from dedicated field campaigns and long-term 
surface networks (Yu et al., 2006 and references therein) has provided strong constraints 
on aerosol optical properties used in satellite aerosol retrievals. Examples of passive 
sensors with the advanced capabilities include the Polarization and Directionality of the 
Earth’s Reflectance (POLDER), MODerate resolution Imaging Spectroradiometer 
(MODIS), Multi-angle Imaging SpectroRadiometer (MISR), and Ozone Measurement 
Instrument (OMI). Major characteristics of these advanced sensors are described in 
Supplementary Online Material (SOM). These sensors provide global measurements of 
AOD over both land and ocean with accuracy on order of 20% or ±0.05 (e.g., Chu et al., 
2002; Remer et al., 2002, 2005; Levy et al., 2010; Hsu et al., 2004; Lyapustin et al., 2011; 
Kahn et al., 2005, 2010; Abdou et al., 2005; Fan et al., 2008; Ahn et al., 2007). 

Remote sensing of AOD from the above daytime polar-orbiting sensors is complemented 
by measurements from infrared sensors such as the Atmospheric Infrared Sounder 
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(AIRS) and geosynchronous satellites. AIRS observations have also been used to detect 
large particles such as soil-derived dust in both day and night (Pierangelo et ah, 2004; 
DeSouza-Machado et ah, 2006; Peyridieu et ah, 2010; DeSouza-Machado et ah, 2010). 
AOD measurements from the polar-orbiting satellites do not sample frequently enough to 
characterize the wide range of daily AOD variability in different locations and seasons 
(Zhang et al., 2012). NOAA’s GOES has been providing 30 min AOD measurements 
over the U.S. for many years (Knapp et al., 2002; Prados et al., 2007). The SEVIRI 
instrument on board Meteosat Second Generation (MSG) also offers new capabilities for 
monitoring aerosol transport over the Atlantic and Mediterranean at high temporal and 
spatial resolution (Thieuleux et al., 2005). Because these sensors generally have coarse 
spatial resolution and limited wavelength coverage, the AOD products from them have 
never reached the same level of accuracy and quality as aerosol products from the EOS- 
era polar orbiting missions. In the coming era of geostationary remote sensing, new 
instruments such as the Advanced Baseline Imager (ABI) onboard GOES-R (scheduled 
for launch in 2015) will provide AOD measurements with significantly improved 
accuracy (Laszlo et al., 2008). 

Observing aerosols above low-level clouds from satellites has also been explored 
recently. CALIOP can measure backscatter and infer aerosol extinction profiles above 
low-level clouds (Winker et al., 2010) or above-cloud AOD and particle properties based 
on lidar observations of depolarization ratio or color ratio for underlying clouds (Hu et 
al., 2007; Chand et al., 2008), albeit with very limited coverage. Several passive sensors 
can also measure AOD above clouds, because the presence of aerosols can significantly 
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modify some attributes of radiances reflected by underlying clouds, such as the polarized 
light in certain scattering-angle ranges (Waquet et ah, 2009) and the spectral contrast of 
radiation (Hennan J. et ah, 1997; Torres et ah, 2007; Fromm et ah, 2008). These 
constitute a basis for retrieving above-cloud AOD using multiangle polarization 
measurements from POLDER (Waquet et ah, 2009) and the backscattered UV (BUV) 
measurements from OMI (Torres et ah, 2012), respectively. A recent analysis shows that 
OMI absorbing aerosol index (AAI) in cloudy conditions (as determined from MODIS 
cloud measurements) correlates well with CALIOP above-cloud AOD, having a 
relationship that depends on cloud albedo (Yu et al., 2012b). The results may suggest the 
potential of combining OMI AAI and MODIS cloud masking to empirically derive 
above-cloud AOD with a spatial coverage much more extensive than CALIPSO (Yu et 
al., 2012b). These above-cloud aerosol measurement capabilities could provide new 
insights into aerosol intercontinental transport, because the transport often occurs above 
clouds (e.g., Dirksen et al., 2009). The above-cloud aerosol observations could reveal to 
some extent the role of cloud processes in pumping atmospheric boundary layer (ABL) 
aerosol up to the free troposphere (FT) for the subsequent cross-ocean transport. They are 
also essential to assessing aerosol direct radiative forcing (Abel et al., 2005; Chand et al., 
2009) and influence on the fonnation and evolution of clouds (Johnson et al., 2004; 
Wilcox, 2010). 

2.2. Emerging climatology of aerosol vertical distribution 

A climatology of aerosol vertical distribution is emerging in recent years due to the 
deployment of space-based lidar technology, in particular the launch of CALIPSO in 
2006 (Winker et al., 2009). Since then, CALIPSO has been collecting an almost 


13 



Yu et al., Satellite Perspective of Aerosol Intercontinental Transport 


continuous record of high-resolution (333 m in the horizontal and as fine as 30 m in the 
vertical), day and night profiles of aerosols and clouds attenuated backscatter at 532 and 
1064 mn wavelengths, along with polarized backscatter in the visible channel, covering 
the globe from 82°N to 82°S (Winker et al., 2009; 2010). Aerosol profiling from 
CALIOP is not limited to cloud-free conditions. Lasers can penetrate thin cirrus clouds 
and sense aerosol layers abvove low-level clouds. As such, a climatology of global 
aerosol and cloud vertical distributions over 5+ years is emerging. 

Passive sensors also show some promise in sensing the height of aerosol layers. 

Radiances at the atmospheric windows (e.g., 8-12 um) are mostly sensitive to the altitude 
of dust layers, due to relatively large particle size and frequently high AOD, compared to 
other aerosol types (Pierangelo et al., 2004; DeSouza-Machado et al., 2006). Having 
nearly complete coverage of the 9-12 pm spectrum with high spectral resolution, AIRS 
has shown capability of detecting dust layer height both day and night (Peyridieu et al., 
2010; DeSouza-Machado et al., 2010). The O 2 -O 2 absorption feature at 477 mn from 
OMI also contains useful infonnation on the altitude of cloud or aerosol layers (Acarreta 
et al., 2004). Scattering by clouds or aerosols reduces the probability of absorption by O 2 - 
O 2 below these layers. The higher the cloud or aerosol layer, the larger the O 2 -O 2 
absorption reduction. From stereo imaging, the heights of clouds and aerosol plumes can 
also be derived geometrically (Moroney et al., 2002; Muller et al., 2002). This technique 
is not dependent on radiometric calibration accuracy, but is effective only near aerosol 
sources, where plumes exhibit sufficient spatial contrast to be matched in multiple 
angular views (Kahn et al., 2007; Scollo et al., 2010; Val Martin et al., 2010). It has been 
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shown that these passive techniques can be used to estimate the altitude of dust and 
smoke plumes under cloud-free conditions, providing essential information that is 
complementary to CALIPSO data (Kahn et ah, 2008; Dirksen et al., 2009; Peyridieu et 
ah, 2010). 

2.3. Enhanced capabilities for characterizing aerosol properties 

Satellite sensors, though not able to directly measure aerosol composition, can constrain 
particle size (i.e., fine versus coarse), shape (i.e., spherical versus non-spherical), and 
absorption (i.e., scattering versus absorbing), from multi- wavelength, multi-angle, and 
polarization measurements. The MODIS multi- wavelength measurements over wide 
spectral range make it feasible to retrieve quantitative aerosol size parameters (e.g., 
effective radius, fine-mode fraction or FMF of aerosol optical depth) over dark water 
(Kaufman et al., 2002; Remer et al., 2005; Levy et al., 2010). Figure 2 shows spatial 
variations of aerosol type on April 11, 2001 as revealed by the composite of aerosol 
optical depth and fine-mode fraction from MODIS observations and the Goddard 
Chemistry Aerosol Radiative Transport (GOCART) simulations. Industrial pollution and 
biomass burning aerosols are predominately small particles and shown as red, whereas 
mineral dust and sea salt are primarily large particles and shown as green. Methods have 
been developed to estimate the optical depth of combustion aerosol (from industrial 
activities and biomass burning) and dust using the FMF and AOD (e.g., Kaufman et al., 
2005a; Yu et al., 2009b). These methods are dependent on several important assumptions 
based on regional climatology and are not the same as using the infonnation content of 
the observations directly to determine aerosol type. However, the methods have provided 
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several important estimates of dust and pollution aerosol transport across ocean basins. 
The use of MODIS deep blue channel also provides an approach of retrieving aerosol 
single-scattering albedo over land (Hsu et ah, 2004). 

The MISR multi-angle data sample scattering angles ranging from about 60° to 160° in 
mid-latitudes, yielding constraints on particle size, single-scattering albedo (SSA) (Kahn 
et al., 1998; 2001; 2005; 2010; Chen et al., 2008) and shape (Kalashnikova and Kahn, 
2006). Having the ability to separate three-to-five bins in aerosol size (e.g., “small,” 
“medium,” and “large”), two-to-four bins in SSA, and spherical vs. non-spherical 
particles under favorable retrieval conditions, MISR can map aerosol-air-mass-types. 
These are derivations directly from the infonnation content of MISR’s observations, 
without the additional burden of the assumptions needed by MODIS to separate particle 
types. Figure 3 shows such two examples that illustrate MISR’s capability of 
characterizing the load and optical properties (size and shape) of smoke and dust plumes. 

POLDER’s unique capability of measuring spectral and angular polarization can yield 
additional insights into aerosol type. When the geometrical conditions are optimal 
(scattering angle ranging between 90°-160°), POLDER can derive the shape of coarse- 
mode particles over ocean (Herman et al., 2005). This allows for a discrimination of non- 
spherical, large dust particles from spherical, large sea-salt particles over ocean. Ongoing 
development of POLDER enhanced algorithms is expected to produce measurements of 
aerosol absorption and altitude (Tanre et al., 2011; Dubovik et al., 2011). OMI AAI can 
be used to distinguish UV-absorbing aerosols such as dust, smoke, and volcanic ash from 
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scattering aerosol like sulfate (Hsu et al., 1996; Herman J. et al, 1997). CALIOP allows 
for the measurements of particle depolarization ratio at 532 mn, a ratio of perpendicular 
component to parallel component of backscatter by aerosol particles, which can be used 
to distinguish non-spherical aerosol (e.g., dust, volcanic ash) from spherical aerosol (e.g., 
industrial pollution, biomass burning smoke, and marine aerosol) (Winker and Osborn, 
1992; Cattrall et al., 2005). These quantities are of interest in-and-of themselves for 
identifying aerosol air mass types, and should also help further refine the accuracy of 
space-based AOD retrievals and particle property determinations. 

The measurements of aerosol microphysical and optical properties can be used in 
conjunction with measurements of trace gases and identified source locations to improve 
aerosol characterization. SCfi and NCfi are important aerosol precursors that are now 
measured by OMI (Carn et al., 2007; Krotkov et al., 2008; Bucsela et al., 2006), the 
Global Ozone Monitoring Experiment (GOME) (Martin et al., 2002; Khokhar et al., 
2005), and the Scanning Imaging Absorption Spectrometer for Atmospheric 
Chartography (SCIAMACHY) (Lee et al., 2008; Richter et al., 2005; Martin et al., 2006). 
The spatial and temporal correlations between concurrent measurements of AOD and 
tropospheric columns of NO 2 and SO 2 provide insights into regional pollution control of 
combustion processes (Veefkind et al., 2011). CO and fine particles have different 
atmospheric lifetimes, but have common sources from industrial pollution and biomass 
burning. Correlations between CO concentration and fine-mode AOD on different spatial 
and temporal scales can provide insights into aerosol emissions, transport, and chemical 
processes (Edwards et al., 2004; Bian et al., 2010). CO is now being derived from the 
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observations of the Measurements of Pollution in the Troposphere (MOPITT) (Deeter et 
ah, 2003), AIRS (Warner et ah, 2007), the Tropospheric Emission Spectrometer (TES) 
(Rinsland et ah, 2006), and the Microwave Limb Sounding (MLS, for upper troposphere 
and low stratosphere (Livesey et ah, 2008). In addition, satellite measurements offer the 
capability of identifying dust and smoke source locations. For example, AI has been used 
to identify major dust source regions (Prospero et al., 2002). Current satellites can detect 
active fires globally on a daily basis (see the review by Ichoku et al., 2012). These source 
identifications should provide additional constraints on microphysically based aerosol 
characterizations . 

2.4. Feasibility of synergistic aerosol characterization 

Rarefy can a single sensor or variable provide unambiguous and complete aerosol 
characterization, because aerosol concentration and composition vary over a wide span of 
temporal and spatial scales. An approach that uses complementary information from 
multiple sensors can lead to improved aerosol characterization. Such synergy is feasible 
because of the recently formed A-Train, a constellation of satellites that overpass the 
equator within a few minutes of each other (Stephens et al., 2002). Radiance 
measurements from multiple A-Train sensors may be combined to retrieve new 
information of particle properties and/or improve retrieval accuracy (e.g., Jeong and Hsu, 
2008; Satheesh et al., 2009). The process is challenging, and a great deal of effort is 
needed to explore the possibilities. 

Currently the most often used multi-sensor synergy is to integrate retrieved products from 
several sensors to get more complete AOD coverage or constraints on aerosol properties 
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(e.g., Anderson et al., 2005; Christopher and Jones, 2007; Kalashnikova and Kahn, 2008; 
Satheesh et al., 2009; Yu et al., 2012b, 2012c). Such a synergy is particularly helpful in 
regions having complex mixtures of aerosol components. For example, the tropical 
Atlantic Ocean is often influenced in boreal winter by a complex mixture of dust and 
smoke transported from North Africa (e.g., Haywood et al., 2008; Tesche et al., 2009), 
which poses a big challenge to aerosol characterization and assessment of aerosol impacts 
(Yu et al., 2009b; Ben- Ami et al., 2012). Figure 4a shows a synergistic characterization 
of dust and smoke over the Gulf of Guinea on January 31, 2008 using a range of 
measurements from multiple A-Train sensors. OMI Al, a good indicator of elevated dust 
and smoke, shows two aerosol plumes, one in the northwestern part of the Gulf and 
another due south toward the coast. For convenience we refer to them as the 
“northwestern regime” and the “southern regime”, respectively. The northwestern regime 
contains aerosol plumes with MODIS AOD of 0.7~0.9, AIRS CO of 2.9~3.1xl0 18 
molecules cm' 2 , and MODIS FMF of up to 0.8. The southern regime is characterized by 
higher AOD, up to 1.5, but relatively low CO (~2.8x l() ls molecules cm 2 ) and FMF 
(0.2~0.6). The FMF seems to anti-correlate with the incremental AOD/CO ratio 
(AAOD/ACO, which is calculated with background CO = 2.4 x 1 0 1 8 molecules cm 2 and 
AOD = 0.1). AAOD/ACO in the southern regime is 2~4xl0" 18 cm 2 molecules' 1 , which is a 
factor of 2-4 larger than that in the northwestern regime. A higher AAOD/ACO indicates 
that not only biomass burning smoke contributes to the AOD enhancement (which would 
cause both AOD and CO to increase), but that dust aerosol also contributes (AOD 
increases but CO doesn’t). These observations seem to consistently suggest that the 
northwestern regime is likely influenced more by fine-mode smoke than dust aerosol, 
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whereas the southern regime is influenced by dust and smoke mixed with a substantial 
dust contribution. The CALIPSO lidar measurements, albeit limited by their cross orbit 
coverage, provide layer structures of the plumes. The plume top in the northwestern 
regime is at about 1 km, much lower than 2-2.5 km in the southern regime. Also, the 
particulate depolarization ratio in the southern regime is somewhat higher than that in the 
northwestern regime, again indicating larger contribution of non-spherical dust in the 
southern regime. These observed distinct plumes over the Gulf can be further traced to 
their source regions through a HYSPLIT backward trajectory analysis ( Draxler and 
Rolph, 2003), as shown in Figure 4b for two receptor locations, i.e., [9°N, 15°W] in the 
northwestern regime and [4°N, 7.5°W] in the southern regime. The air mass arriving at 
both locations on January 31 , 2008 1400UTC had passed over the smoke-active Sahel 
region during the January 28-30 period (Tesche et al., 2009). During the January 26-27 
period, the air mass arriving at [4°N, 7.5°W] passed over center of the Bodele Depression 
- the major source of Saharan dust in winter (Koren et al., 2006) and may have picked up 
a significant amount of dust, whereas that arriving at [9°N, 15°W] drifted north of the 
Bodele Depression. This case study demonstrates the advantage of multi-sensor synergy 
in characterizing aerosol type and long-range transport. Column-integrated aerosol 
properties observed by passive sensors are complemented by vertical layer structures 
observed by the active sensor. Aerosol type characterization based on particle size and 
shape is cross-examined with CO measurements and traced to source locations through 
the backward trajectory analysis. 
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The combination of polar-orbiting and geostationary satellite data could provide sub- 
daily infonnation on the evolution of aerosol plumes during transport. Polar-orbiting 
satellites can sample a given location no more than twice a day, which is insufficient to 
capture diurnal or daytime variations adequately (Zhang et ah, 2012). A use of 
geostationary satellite observations, such as half-hourly GOES AOD measurements, can 
characterize the daytime variations of intercontinental transport events (Kondragunta et 
al., 2008; Laszlo et al., 2008). Geostationary satellites can also capture diurnal variations 
of biomass burning smoke emissions with hourly or sub-hourly frequency (Zhang and 
Kondragunta, 2008; Reid et al., 2004), essential inputs to chemical transport models. 

In summary, major advances have been achieved in both passive and active aerosol 
remote sensing from space in the last decade or so, providing better coverage, higher 
spatial resolution, improved retrieved AOD accuracy, and constraints on particle 
properties. Space-based lidars are for the first time providing global constraints on 
aerosol vertical distribution, and multi-angle imaging is supplementing this with maps of 
plume injection height in aerosol source regions. Satellite remote sensing of aerosols 
above clouds has also been demonstrated. Major advances have been made in 
distinguishing aerosol types from space, and the data are now useful for evaluating 
aerosol transport model simulations of aerosol air mass type distributions and transports, 
particularly over dark water. These advances provide unprecedented opportunities to 
develop a satellite-based approach, supplemented by in situ measurements and model 
simulations, to assess aerosol intercontinental transport. In the next section, we 
summarize major efforts and findings over the last decade and discuss future efforts at 
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maximizing the use of satellite observations to assess aerosol ICT and impacts on climate 
and air quality. 

3. Application of Satellite Measurements to Characterizing and Estimating Aerosol 
Intercontinental Transport 

The accumulation of satellite observations over the past decades has significantly 
advanced the understanding of aerosol intercontinental transport and its impacts on global 
and regional climate change and air quality. In this section we review how these 
observations have been employed to characterize the details of aerosol ICT episodes 
(3.1), study seasonal and inter-annual variations of aerosol ICT (3.2), and estimate the 
material fluxes of aerosols leaving a continent and arriving at another (3.3). In addition to 
these standalone observation-based applications, satellite measurements have also played 
an important role in evaluating and constraining model simulations of aerosol 
intercontinental transport (3.4). Both progress and challenges are discussed. 

3.1 . Episode-based characterization of aerosol ICT from satellites 

The daily or hourly observing frequency of satellites allows for tracking the evolution of 
significant aerosol ICT episodes, such as dust storm outbreaks, severe forest fire smoke, 
and abrupt volcanic eruption plumes. The extensive spatial coverage of satellite 
observations also helps establish linkages between in situ measurements in different 
locations that might be thousands of kilometers apart (Ansmann et al., 2009; Ben- Ami et 
al., 2010). 
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One of the first uses of satellites to characterize aerosol ICT episodes was simply to 
provide imagery of major aerosol events as those plumes crossed oceans. For example, 
Husar et al. (2001) used a combination of SeaWiFS, GOES 9, GOES 10 and TOMS UV 
aerosol index (UVAI) to analyze by hand the approximate location of a dust plume that 
developed over the Asian Gobi desert and then was transported across the Pacific Ocean 
over the course of 9 days. The arrival of the dust in the west shore as determined by 
satellite analysis, was consistent with elevated concentrations of dust particulate matter 
measured at in situ stations in California and Oregon. Although this satellite study and 
other similar ones provided little quantitative infonnation about the aerosol transport, it 
was convincing illustration of aerosol transport across ocean basins. 

Satellite analysis provided a more quantitative characterization of aerosol ICT in the 
1980s and 1990s when applied to major volcanic eruptions that caused sulfate aerosol to 
form in the stratosphere. Solar occultation instruments such as the Stratospheric Aerosol 
and Gas Experiment I and II (SAGE and SAGE II), the Stratospheric Aerosol 
Measurements II (SAM II) or the Cryogenic Limb Array Etalon Spectrometer (CLAES) 
were able to track the spread of the El Chichon and Mount Pinatubo aerosol plumes as 
they spread from the point of injection above the volcano, eventually covering the entire 
tropical belt over the course of several years (Lambert et al., 1997; Bauman et al., 2003). 
These observations provided monthly and zonally averaged measures of aerosol 
extinction in the stratosphere, sometimes resolved into vertical profiles. Likewise the 
Advanced Very High Resolution Radiometer (AVHRR) followed the Mount Pinatubo 
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plume as long as it dominated the total column aerosol loading (Long and Stowe, 1994; 
Stowe et al., 1997; Stowe et al., 2002). 

In contrast the EOS-era instruments offer a higher degree of quantitative characterization 
of episodic transport. For example, making use of the Terra-MODIS sensor’s wide 
swath, validated AOD product and information that proportions AOD into fine and 
coarse modes, daily snap shots of the Earth’s aerosol system became available (Remer et 
al., 2005). Applying temporal averaging to these daily global images and accumulating a 
year’s set of images, animations were constructed that clearly showed intercontinental 
transport of fine and coarse particles [see Supplementary Online Material]. Saharan dust 
crossed the Atlantic not as a continuous stream, but in waves. Similarly, specific 
episodes of fine mode aerosol left the southeast U.S. for Europe and multi-modal aerosol 
crossed the Pacific as well-defined events. 

However, the more significant EOS-era improvement was the ability to characterize 
particle properties, even beyond the simple designation of fine- or coarse-dominated 
AOD. Now satellites had opportunities for not only tracking the spatial extent of the 
influence of episodes, but also of characterizing particle property evolution during the 
cross-ocean transport. Kalashnikova and Kahn (2008) analyze the trans-Atlantic 
transport of four dust plumes during a peak dust season using MISR and MODIS aerosol 
products. They report the evolution of horizontal extent, AOD, particle size, and particle 
shape within the plumes. As dust is transported across the Atlantic, AOD and its non- 
spherical fraction systematically decrease, whereas Angstrom exponent increases. For the 
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thick part of the plume (AOD > 0.5), Angstrom exponent and spherical fraction of AOD 
increase by about 30% as the plumes travelled from the coast of North Africa to the 
Caribbean Sea. The study also demonstrates how a combination of complementary 
information from the two sensors can be used to gain better characterization of cross- 
ocean dust transport. MODIS provides more extensive coverage, whereas MISR’s multi- 
angle retrievals include dust properties and fill in areas where glint precludes MODIS 
retrievals. 

The most prevalent use of EOS-era satellite sensors to characterize aerosol ICT events or 
episodes makes use of the CALIPSO lidar that offers vertical profiling for identification 
of transport height, color ratio that provides qualitative particle size and depolarization 
that identifies non-spherical particles such as dust. With successive CALIPSO lidar 
profiles, Liu Z. et al. (2008) track the cross-Atlantic transport of a dust plume that 
originated in the Sahara desert on August 17, 2006 for a period of 10 days. They found 
that it took about a week for the dust plume to reach the Caribbean from North Africa, as 
the optical depth decreased significantly and the plume top descended by about 3.3 km. 
However, unlike the plumes studied with MISR observations (Kalashnikova and Kahn, 
2008), the size and shape of the dust particles didn’t show discernable changes during 
transport. After entering the Caribbean region, the transport slowed, and dust was mixed 
with other types of aerosol over the Gulf of Mexico, as indicated by changes in color 
ratio and lidar ratio. 
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Trans-Pacific transport of dust has also been characterized by using successive CALIOP 
curtains in conjunction with other satellite sensors or models, including Huang et ah, 
(2008), Uno et al. (2008), Eguchi et al. (2009), Uno et al. (2009), Yumimoto et al. (2009), 
and Uno et al. (2011). Figure 5 shows a three-dimensional view of the hemispherical 
journey of strong dust event originated from the Taklimakan desert during May 09-22, 
2007 as provided by CALIOP for vertical distributions of depolarization ratio (a) and 
extinction coefficient (b), and a global aerosol transport model for horizontal distribution 
of dust extinction coefficient (c) (Uno et al., 2009). The HYSPLIT trajectory also 
illustrates the dust transport path in both horizontal and vertical directions. For this 
extreme dust storm, dust was lofted to 8-10 km and was transported to the mid-Pacific in 
4 days, and in about 13 days, it traveled more than one full circuit around the Northern 
Hemisphere. On returning to the northwestern Pacific Ocean, the dust then descended 
within a large-scale high-pressure system. 

Strong smoke and volcanic eruption events have also been characterized by satellite 
measurements. Dirksen et al. (2009) use OMI and CALIOP measurements to track a 
smoke aerosol plume released by intense forest fires in Australia that circumnavigated the 
world in 12 days. The study highlights the need to include the effect of pyro-convective 
lofting and the effect of smoke absorption on the atmospheric thennodynamics in global 
models. CALIPSO lidar data have also been used to characterize the transport of ash 
from the Eyjafjallajokull volcano in April 2010, including height and vertical extent of 
the layers and particle shape during both day and night (Winker et al., 2012). The 
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observed particle layers have strong lidar depolarization similar to mineral dust, 
suggesting the predominance of ash over sulfate in the plumes (Winker et ah, 2012). 

The transport of aerosols from mid-latitudes to the polar regions has great implications 
for climate change. Satellite measurements, in particular those from the active sensor 
CALIOP over day and night, are useful to some extent for characterizing transport events 
of eastern Asian aerosols to the Arctic (Di Pierro et al., 2011) and of Patagonian dust to 
Antarctica (Gasso and Stein, 2007; Li et al., 2010). However, the usually low aerosol 
loading in polar regions requires a space-borne lidar with a detection limit lower than 
CALIOP. In comparison to active sensors, passive sensors have much greater difficulties 
in detecting the poleward transport of aerosols at high latitudes, because of low solar 
signal and hence more uncertain aerosol retrieval, and the presence of abundant clouds in 
the roaring forties (Gasso and Stein, 2007; Li et al., 2010). 

3.2 Climatological characterization of aerosol ICT from satellites 

The long-term duration of satellite measurements makes it feasible to derive the seasonal 
and inter-annual variations of aerosol ICT. Satellite aerosol products are easily 
accumulated and averaged, and when these statistics are derived over oceans, we assume 
that the elevated aerosol loading is due to transport from continents. However, without 
additional infonnation, either from tracking specific events, analysis of ancillary 
meteorological data, or model trajectories, the statistics of oceanic aerosol retrievals 
cannot be directly attributed to transport from a specific location. It helps if satellites 
provide additional infonnation about the particle properties to identify certain mean 
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aerosol loading to be associated with an aerosol type originating from continents (e.g., 
dust or combustion particles.) 

Seasonal variations of trans-Atlantic dust transport have been extensively studied using 
measurements from passive sensors, often supplemented with ancillary data, (e.g., 
Kaufman et al., 2005c; Schepanski et al., 2009; Huang et al., 2010; Ben-Ami et al., 2012; 
Ridley et al., 2012) and active sensors (Generoso et al., 2008; Ben-Ami et al., 2009; 
Huang et al., 2010; Yu et al., 2010; Ridley et al., 2012). An analysis of aerosol 
observations from TOMS, AVHRR, and MODIS looked for intra-seasonal variations in 
AOD over the tropical Pacific and Indian oceans, with inconclusive results due to 
difficulties in interpreting the satellite products on relatively short spatial and temporal 
scales in these regions (Tian et al., 2008). However, a follow-up study applied the same 
technique over the tropical Atlantic, and showed that about 25% of the total variance of 
MODIS AOD can be accounted for by the intra-seasonal variability related to the 
Madden- Julian Oscillation (MJO) (Tian et al. 2011). Given that dust and smoke are often 
mixed in winter and spring (Haywood et al., 2008), several studies have attempted to 
separate dust from other aerosols in their analysis (Kaufman et al., 2005c; Liu D. et al., 
2008; Huang et al., 2010; Yu et al., 2010; Ben-Ami et al., 2012). This separation can be 
done by using depolarization from CALIOP (Liu D. et al., 2008), by retrieving detailed 
aerosol particle properties from multi- angle and/or polarization measurements (Kahn et 
al., 2007; Hennan et al., 2005) or by combining size parameters derived from satellite 
sensor with assumptions gleaned from climatological analysis (Kaufman et al., 2005a, 
2005c; Yu et al., 2009b). 
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The lidar depolarization method makes use of the fact that spherical particles such as 
pollution, smoke and sea salt do not create a depolarization signal, but irregularly shaped, 
large dust particles do (McNeil and Carswell, 1975). It has been observed that mineral 
dust has a high depolarization ratio of 0. 2-0.4 (Bamaba and Gobbi, 2001; Murayama et 
al., 2001; Liu et al., 2002; Maths et al., 2002; Liu Z. et al., 2008), whereas other types of 
aerosol have near zero depolarization ratios. Therefore this depolarization signal allows 
the lidar to determine dust layers in the vertical profile (Winker and Osborn, 1992). 

Dust can also be separated from combustion particles such as smoke and pollution when 
a sensor characterizes the multi-dimensional reflectance field, including 
multiplewavelengths across a broad spectral range, multiple angle over a wide range of 
scattering angles and polarization states (Dubovik and King, 2000; Mishchenko et al., 
2005). Different subsets of the full characterization offer different capabilities. MISR 
provides multi-wavelength, multi-angle measurements and provides retrievals of particle 
properties of which nonsphericity is the most important for identifying dust 
(Kalashnikova and Kahn, 2006; Kahn et al., 2007). PARASOL-POLDER is also multi- 
wavelength, multi-angle instrument and adds polarization, making it also capable of 
directly retrieving sufficient information on particle properties to discern dust (Herman et 
al., 2005). 

On the other hand, MODIS with its single angle view of each scene and no polarization 
does not measure sufficient information to directly retrieve the particle properties that 


29 



Yu et al., Satellite Perspective of Aerosol Intercontinental Transport 

unambiguously identify dust. Following Kaufman et al. (2005a), the MODIS studies 
assume: 



where Eq. 1 assumes that the total AOD (x) retrieved by MODIS is composed of a dust 
component ( Xd u ), a combustion aerosol component ( x co ) and a marine aerosol 
component ( x ma ). Eq. 2 assumes that the fine mode AOD (Tf me = f x, the fraction of the 
total AOD contributed by fine particles) is composed of a sum of the fine mode AODs of 
each component in Eq. 1 . Note that the fine mode fraction (f) is the ratio of fine mode 
AOD and the total AOD, 

/ = ^ (3) 

T 


MODIS derives two pieces of information, the total AOD (x) and the fine-mode fraction 
(f). This means that all but two of the parameters on the right hand side of Eqs. (1) and 
(2) must be assumed. The assumed quantities are fd u , f co , La and x ma . The values for 
these assumptions are taken from MODIS retrievals of f and x in regions that are 
expected to be pure dust, pure combustion aerosol etc. f ma was found to vary spatially and 
seasonally (Yu et al., 2009b) and is therefore a function of location and time. Given 
these assumptions and the MODIS direct retrievals of f and x, Xd u and x co can be inferred. 
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Note that at no point in this method is it assumed that dust is pure coarse mode (fd u ^ 0) 
or combustion aerosol is pure fine mode (f co ^ 1), but the actual values are fixed for a 
given season and region, and are obtained as described above. 

Because each sensor, to differing degrees of quantification accuracy, can discern dust 
from the total aerosol signal observed, all have been used to characterize seasonal and 
inter-annual variations of dust transport. As an illustration, Figure 6 shows a multi- 
sensor (including MODIS, MISR, and PARASOL) view of seasonal variations of trans- 
Atlantic dust transport and its comparison with the Goddard Chemistry Aerosol Radiation 
Transport (GOCART) model. In general, these satellite observations and the model show 
consistent seasonal variations of trans-Atlantic dust transport, i.e., stronger in summer 
and spring than in winter and fall. On the other hand, notable differences exist in the 
magnitude of dust AOD. For example, although PARASOL only senses the coarse-mode 
non-spherical particles in this plot (Tanre et al., 2011), its AOD value is comparable to 
the sum of fine-mode and coarse-mode dust AOD as seen by MODIS and MISR. Note 
that fine-mode dust makes significant (e.g., 37%) contributions to the total dust AOD for 
MODIS observations (Yu et al., 2009b). In the Gulf of Guinea, MODIS and PARASOL 
show much higher dust AOD than MISR and GOCART do in winter. The complexity of 
dust and smoke mixture in the region creates a challenging environment for 
distinguishing dust from smoke using satellites (Kahn et al., 2009; Yu et al., 2009b). In 
such a situation, MISR would be the sensor of choice, though the known issues with that 
product (e.g., Kahn et al., 2010) should also be taken into account. As PARASOL 
algorithms are further developed to make better use of the full complement of its 
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measurements and information content (Dubovik et al., 2011), we anticipate that its 
retrievals will draw closer to MISR’s. 

In comparison to the satellite observations, GOCART modeled cross-ocean transport 
generally drifts southward in summer and decays more rapidly (i.e., weak transport 
efficiency). Such differences are apparently not limited to the GOCART model, as 
several other models show similar deviations from satellite observations (e.g., Generoso 
et al., 2008; Huneeus et al., 2012; Ridley et al., 2012). Reasons for such observation- 
model discrepancies are not yet fully understood. It has been argued that model treatment 
of aerosol removal processes needs to be improved (Generoso et al. 2008; Prospero et al., 
2010; Huneeus et al., 2012; Ridley et al., 2012). Dust can influence the cross-ocean 
transport, through its feedbacks on atmospheric circulations and the water cycle, that are 
not accounted for in most CTMs. For example, Lau et al. (2009) found that dust 
absorption can enhance rainfall and cloudiness over the West Africa and Eastern Atlantic 
Intercontinental Convergence Zone (ITCZ), but suppresses rainfall and cloudiness and 
rainfall over the West Atlantic and Caribbean region. Wilcox et al. (2010) found that 
Saharan dust outbreaks coincide with a northward shift of the North Atlantic Ocean ITCZ 
in summer. If these feedbacks and coincidences are supported by independent 
observations and can be reproduced by aerosol models, the modeled summertime cross- 
ocean transport would have shifted northward and the transport efficiency or the gradient 
of dust AOD between the East and West Atlantic would increase. 

The seasonal variations of columnar AOD observed by passive sensors are 
complemented by lidar measurements of aerosol vertical distributions. Recently available 
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CALIPSO data confirm the previously observed seasonal differences of trans-Atlantic 
dust transport heights, namely higher in summer than in winter (Generoso et ah, 2008; 
Ben-Ami et ah, 2009; Huang et ah, 2010; Yu et ah, 2010; Ridley et ah, 2012). These 
same analyses also show that a significant part of the dust is transported near and within 
the marine boundary layer, suggesting significant dust deposition to the ocean and hence 
important implications for ocean biogeochemistry. Vertical structures of dust and 
pollution transport on global and hemispherical scales have also been analyzed using 
CALIOP measurements (Liu D. et al., 2008; Yu et al., 2010; Ford and Heald, 2012; Koffi 
et al., 2012). 

Decadal-scale satellite measurements of aerosols, many relying on pre-EOS-era satellites 
to obtain longer time series, have been used to study inter-annual variations of aerosol 
ICT and their controlling factors. Such studies are not necessarily trend studies, which 
are particularly sensitive to calibration drift in the sensors (Levy et al., 2010; Zhang and 
Reid, 2010). If focused on a particular region with high aerosol loading and strong 
interannual amplitudes, the aerosol variability can be distinguished above sensor artifacts 
(Zhang and Reid, 2010). Several studies have investigated interannual variation in the 
export of Saharan dust using satellite data. Moulin et al. (1997) used daily Meteosat 
observations of dust from 1984-1994 and found that inter-annual variability of dust 
export to the North Atlantic Ocean and the Mediterranean Sea are well correlated with 
the North Atlantic Oscillation (NAO). By using multi-decadal TOMS Al and ground- 
based dust measurements, Chiapello et al. (2005) found that the Sahel drought has a 
large-scale impact on dust emissions and trans- Atlantic transport in both the subsequent 
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winter and summer, and that the NAO controls the northern branch of wintertime 
transport. Doherty et al. (2008), also using the multi-decadal TOMS AI data, found that 
African dust is transported into the Caribbean via two routes associated with distinct 
source regions and controlled by different meteorology. In the northern mode, dust 
mobilized from the Sahara is transported westward, controlled primarily by the 
longitudinal displacement of the Azores High. In contrast, in the southern mode, dust 
originating from the Sahel region is transported to the Caribbean via the Gulf of Guinea, 
controlled primarily by longitudinal displacement and surface pressure fluctuations of the 
Hawaiian High. Note that the pre-EOS era sensors of Meteosat and TOMS sensors alone 
cannot distinguish dust from smoke in the above studies, whereas it is necessary to assess 
how inter-annual variations of smoke aerosol may have contributed to the analyzed inter- 
annual variations of dust. 

Climatology of aerosol parameters has been constructed from the EOS-era satellites 
(Remer et al, 2008; Dey and Di Girolamo, 2010 ; Tanre et al., 2011). These do show 
interannual variability of the aerosol system, including interannual variability of aerosol 
distribution across ocean basins. Nevertheless the question of the variability of ICT 
specifically is not adequately addressed in these relatively short time series. 

Undoubtedly, an aerosol data record longer than the 10 years of EOS-era satellite 
observations is essential for a better understanding of the climate impacts of aerosol ICT 
and controlling processes. There have been a few attempts at extending the relatively 
short EOS era datasets by using pre-EOS sensors as proxies. Evan and Mukhopadhyay 
(2010) attempted to extend the MODIS-estimated dust AOD (Kaufman et al., 2005c) 
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over the northern tropical Atlantic back to 1955 by simply scaling AVHRR total AOD 
with the MODIS dust AOD fraction and using a proxy record of atmospheric dust 
(namely crustal helium-4 flux from a coral 5 m under the water near Sal Island). 

However, such an extension would be more credible if build upon understanding and 
reconciliation of AOD differences between the EOS-era and historic sensors (Jeong et ah, 
2005). For this a closer examination of relevant issues associated with individual sensors 
is urgently needed, including instrument calibration, algorithm assumptions, cloud 
screening, data sampling and aggregation, among others (Zhao et al., 2008). Substantial 
efforts have been dedicated toward addressing these issues across the EOS-era sensors 
(Zhang et al., 2005; Kahn et al., 2009; Levy et al., 2009; Zhang and Reid., 2010). 
However, the cross examination of EOS-era sensors and historical sensors remains 
basically unexplored. A reliable long-term data record will only be accessible after these 
issues have been adequately addressed. 

3.3. Estimates of aerosol export and import mass flux from satellites 

The vast majority of studies that address aerosol ICT using satellite data alone have fallen 
under the subject of the previous two sections: characterizing specific events and 
episodes, or providing insight into seasonal and interannual variation. These studies 
progressed from qualitatively following an aerosol plume across ocean basins, to 
gradually adding quantitative information about the evolution of particle properties and 
height. The aerosol climatology across ocean basins could be interpreted as aerosol 
transport, with the relative changes from season-to-season or year-to-year identified, even 
if the quantitative total transport was not addressed. In recent years, a few studies have 
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taken advantage of the advances in aerosol remote sensing to quantitatively estimate the 
trans-Atlantic transport of dust mass (Kaufman et ah, 2005c; Koren et ah, 2006), trans- 
pacific transport of combustion aerosol mass (including industrial pollution and biomass 
burning smoke) (Yu et al., 2008) and dust mass (Yu et al., 2012c), and cross- 
Mediterranean Sea transport of European pollution mass (Rudich et al., 2008). 

Such estimates of material fluxes in continental outflows and inflows generally follow a 
three-step approach. First, satellite measurements of AOD and aerosol microphysical 
properties (e.g., size, shape) are used to distinguish mineral dust from combustion-related 
aerosol (Kaufman et al., 2005a; Yu et al., 2009b; and Section 3.2). Second, AOD for dust 
or combustion aerosol is converted to aerosol (dry) mass concentration profile using in- 
situ measurements of aerosol hygroscopic property (depending on relative humidity and 
aerosol type) and lidar measurements of aerosol vertical distributions (Yu et al., 2008, 
2012). Third, the aerosol mass concentration in combination with zonal wind speed is 
used estimate the zonal export and import of aerosol. 

There has been a progression in these studies as new infonnation became available. The 
first study of this type (Kaufman et al., 2005c) determined aerosol height from a few field 
experiments, confirmed by applying a technique that used paired Terra and Aqua MODIS 
images to determine transport displacement over a few hours and correlated this with 
wind speed at various heights using reanalysis data (Koren and Kaufman, 2004). They 
also assumed single value for f ma . fd u , and f co . In contrast, the most recent study of this 
type (Yu et al., 2012c) uses CALIPSO lidar-derived climatology to determine transport 
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height, and accepts that f ma varies as a function of latitude and season (Yu et al. 2009b). 
Furthermore, the backbone of these studies is the MODIS aerosol AOD, which has 
undergone considerable refinement over the 7-year span of these studies, including the 
option of selecting a severely cloud-cleared and uncontaminated product (Shi et ah, 2011) 
that was used by Yu et al. (2012c) to correct cloud contamination in AOD. 

The satellite-based estimates of aerosol import and export in the zonal direction have 
offered important insights into aerosol intercontinental transport and its impacts on 
regional air quality, climate change, and biogeochemical cycle. Kaufman et al. (2005c) 
estimated that out of the 240±80 Tg of dust transported annually from Africa to the 
Atlantic ocean, 140±40 Tg is deposited in the Atlantic ocean and 50± 15 Tg reaches and 
fertilizes the Amazon basin. This estimated dust input to the Amazon basin is four times 
the previous estimate by Swap et el. (1992), which might explain a paradox regarding the 
source of nutrients for the Amazon forest. Koren et al. (2006) further estimated that about 
half the annual dust supply to the Amazon basin is emitted from a single, small-area spot, 
namely the Bodele depression, located northeast of Lake Chad. Recently, model 
sensitivity studies have argued that these estimates may represent the upper bound 
(Schepanski et al., 2009; Ridley et al., 2012), because of uncertainties associated with the 
assumed aerosol vertical distributions (Schepanski et al., 2009). On the other hand, model 
simulations tend to give a lower bound, due to excessive wet removal in models (Ridley 
et al., 2012; Huneeus et al., 2012). 
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Trans-Pacific aerosol transport and its contribution to North America have been assessed 
based on satellite measurements by Yu et al. (2008) and Yu et al. (2012c). Figure 7 
shows the satellite-based estimate of meridional variations of annual trans-Pacific dust 
and combustion aerosol import to North America (at the west shore) and its comparisons 
with GOCART and the Global Modeling Initiative (GMI) model simulations. The 
satellite-based estimate of trans-Pacific import is 60 Tg yr' 1 with an estimated uncertainty 
of 55-100%, which is predominated by dust (about 90%) entering the west shore between 
30°-50°N. The GOCART and GMI model is respectively lower than satellite-based 
estimate by 15% and 42%. The much lower import by GMI presumably results from very 
efficient dust scavenging by convection in the model. Also the satellite-based estimate is 
generally higher than model simulations in regions south of 50°N, but are lower north of 
50°N. By taking into account the trans-Atlantic dust import (north of 20°N) of 4 Tg yr' 1 
(Kaufman et al., 2005c), the total aerosol import to North America is estimated at 64 Tg 
yr' 1 , which is comparable to the estimated total (69 Tg yr' 1 ) of domestic emissions and 
production of aerosols (Yu et al., 2012c). The air quality impact of imported aerosols 
may be minor, because the transport occurs predominantly above the boundary layer, as 
indicated by CALIPSO lidar measurements (Yu et al., 2012c). However, the imported 
dust and combustion aerosol can have significant impacts on weather and climate. Yu et 
al. (2012c) estimated that dust and combustion aerosols from intercontinental transport 
collectively introduce a reduction of cloud-free net solar radiation of -1.7 and -3.0 Wm" 2 
at top-of-atmosphere (TOA) and surface, respectively, which represents 31% and 37% of 
the total direct aerosol radiative effect over North America. For comparison, a recent 
model simulation estimates that the decrease of US anthropogenic aerosols over 1990- 
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2010 yields a decrease of 0.8 Wm~ 2 for all-sky direct radiative forcing at TOA 
(Leibensperger et al., 2012). Besides the aerosol direct radiative effects, the imported 
aerosols could have other significant impacts, such as changing atmospheric stability by 
absorbing solar radiation (Yu et al., 2002; Ramanathan and Cannichael, 2008), altering 
cloud and precipitation processes through acting as ice nuclei (Sassen et al., 2002; Ault et 
al., 2011), and accelerating the melting of snow in the Sierra Nevada by deposition on 
snow (Hadley et al., 2010). It is possible that a change in meteorological conditions 
resulting from these processes could further influence surface PM concentrations and air 
quality. These potential impacts need to be investigated with the state-of-the-science 
modeling. 

The dust deposition into oceans has important implications for ocean biogeochemistry 
and climate (Jickells et al., 2005; Mahowald et al., 2005). Considering that aerosol 
meridional transport is generally much weaker than zonal transport, a difference between 
zonal export and import could be considered as a proxy for aerosol depositional flux into 
the ocean during the cross-ocean transport. The ratio of import to export flux is generally 
referred to as “transport efficiency”. For the same export in continental outflow, higher 
transport efficiency generally corresponds to lower aerosol deposition over ocean. Table 
2 compares satellite-based estimates of the trans-Pacific dust transport efficiency and dust 
deposition to the North Pacific (Yu et al., 2012c), with results based on model 
simulations and in situ observations in the literature. The annual trans-Pacific dust 
transport efficiency is estimated at about 40% from MODIS and CALIOP measurements 
(Yu et al., 2012c), which is at the upper bound of the model simulations listed in Table 2 


39 



Yu et al., Satellite Perspective of Aerosol Intercontinental Transport 


(15-43%). Note that the highest model value of 43% is estimated for an extreme event, 
where dust was transported rapidly at 8-10 km elevation (Uno et ah, 2009), whereas the 
lowest model value of 15% was estimated for dust transport in summer, when dust 
removal from atmosphere is faster and dust transport is slower than in other seasons. The 
trans-Pacilic dust transport efficiency is also reasonably higher than the 26% for trans- 
pacific pollution aerosol estimated from satellite (Yu et al., 2008), because dust is 
transported at generally higher altitudes than the pollution aerosol (Yu et al., 2012c). 
Satellite-estimated dust deposition to the North Pacific of 84 Tg a' 1 is compared to the 
estimated range of 31-96 Tg a" 1 in the modeling literature. (These values exclude the 
highest and lowest estimates from two early studies, as well as two estimates made only 
over the East China Sea.) 

The satellite-estimated mass fluxes are subject to large uncertainties, on an order of 55- 
100% (Yu et al., 2008; 2012c). Major sources of uncertainty include satellite-derived 
component AOD, mass extinction efficiency, hygroscopic properties, transport height, 
and satellite sampling (Yu et al., 2008; Kaufman et al., 2005c). To accurately assess the 
potentially large impacts of aerosol ICT, the uncertainties represented by the above 
studies must be significantly reduced. The component AOD derived from MODIS 
measurements entails significantly larger uncertainties than the total AOD. As discussed 
in section 2 and demonstrated by Kalashnakova and Kahn (2008), by combining the 
MODIS observations with aerosol type retrievals from other sensors such as CAILPSO, 
AIRS, or MISR, component AOD can be derived with better accuracy. 
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Determining seasonally and geographically dependent aerosol vertical profiles or 
transport heights is another major source of uncertainty in estimating the aerosol mass 
fluxes (Yu et al., 2008; Schepanski et ah, 2009). To calculate mass fluxes, vertical 
profiles are needed for selecting representative wind speed, and representative RH to 
account for particle growth by humidification. The sensitivity of aerosol mass flux to 
transport height varies with region and season, depending on wind shear patterns and 
atmospheric humidity profiles. The CALIPSO lidar measurements have recently been 
used to determine the climatology of aerosol vertical distributions over broad areas with 
seasonal distinctions (e.g., Liu D. et al., 2008; Huang et al., 2010; Yu et al., 2010, 2012c), 
which provide tighter constraints on the estimated material fluxes. Because of the near- 
zero swath of CALIPSO lidar, however, it is difficult to examine daily variations of 
aerosol fluxes. Because of CALIPSO’s low sensitivity to vertically resolved extinction, 
some aerosol plumes having moderate AOD but large vertical extent may be missed by 
the lidar. This drawback can be mediated by suborbital measurements at key locations to 
calibrate the satellite data sets. 

Several satellite-sampling issues also contribute to the uncertainties in the flux estimates. 
The above studies have only used polar-orbiting satellite measurements, having sampling 
frequency of no greater than twice per day. An integration of polar orbiting and 
geostationary satellites could be pursued to characterize the diurnal variations of aerosol 
export and import. Current studies have also assumed that dry aerosol mass derived from 
satellite retrievals in cloud-free conditions are representative of all-sky conditions, even 
over frontal cloud systems. Unlike in broken cloud fields where satellite retrievals of 
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aerosols are sometimes possible in the gaps between clouds, frontal cloud systems 
usually extend over wide areas and completely eliminate aerosol retrievals in those areas 
using current approaches. For example, because the trans-Pacific transport is associated 
with the wann conveyor belt (WCB), and the WCB is associated with rising air in warm 
fronts that also produce widespread cloudiness, passive sensors such as MODIS can miss 
important transport events. If dry aerosol mass in these events is significantly different 
from that observed under cloud-free conditions, the estimated aerosol flux could have 
large uncertainties or biases. It has been suggested that a variety of compensating 
processes in and nearby clouds control the fonnation and removal of aerosols and hence 
detennine their loading. While clouds can effectively remove aerosols through 
scavenging and rainout, they can also generate aerosol (e.g., in-cloud aqueous production 
of sulfate). Some studies suggest that significant biases may exist in some regions, when 
satellite measurements only for cloud-free conditions are assimilated into a model (Reid 
et al., 2004; Zhang and Reid, 2009). Recent developments in satellite retrievals as 
discussed in section 2.2, could provide partial constraints on the aerosol load above 
clouds. 

3.4. Evaluation and constraint of model simulations by satellite measurements 

The previous sections addressed the contribution that satellites make towards 
characterizing aerosol ICT when satellite data is used exclusively. Satellites also make a 
contribution when they work in support of model simulations. The satellite role in 
support of modeling efforts of intercontinental aerosol transport takes generally three 
paths: (1) characterization of the aerosol sources including location, timing, injection 
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height, and magnitude of emissions, (2) evaluation of aerosol transport simulations, and 
(3) assimilation of aerosol products into the model for improving model forecasts. 

3.4.1 Source characterization 

For a model to properly simulate aerosol distribution and transport, it must first 
incorporate accurate location, timing, and magnitude of particles and gaseous precursor 
emissions. In comparison to industrial emissions, episodic emissions such as biomass 
burning smoke from fires and mineral dust from dust storms are more difficult to 
accurately represent in a model. Correct transport modeling of smoke aerosol is almost 
entirely dependent on satellite-derived products that either identify and count fire 
hotspots or more quantitatively measure fire radiative power and relate that to aerosol 
emissions (Ichoku et al., 2008; Pereia et al., 2009; Reid et al., 2009; Giglio et al., 2010; 
van der Werf et al., 2010; Ichoku et al., 2012 and references therein; Petrenko et al., 
2012). The diurnal pattern of smoke emissions can be detennined by applying overpasses 
at multiple times per day by the twin MODIS sensors on the polar orbiting Terra and 
Aqua satellites (Vennote et al., 2009; Ichoku et al., 2008), or using geostationary satellite 
observations (Zhang and Kondragunta, 2008). For mineral dust, most models rely on 
satellite data of land surface classification to identify the location of deserts, and a few 
models use satellite vegetation index data to impose the seasonal variation of the surface 
bareness for better temporal variation of dust emission (e.g., Zender et al., 2003; Kim et 
al., 2012). More direct use of satellite AOD for dust emission has merged, such as the 
AOD distributions from MISR (Koven and Fung, 2008), Meteosat SEVIRI (Schepanski, 
et al., 2007) or more recently MODIS Deep Blue (Ginoux et al., 2012). Satellite 
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observations have also been used to cap the total emissions of carbonaceous aerosol 
(Heald et ah, 2010). 

Satellites also contribute directly to aerosol source characterization through inverse 
modeling techniques. In these techniques an observed satellite field such as satellite- 
derived AOD or measured radiances are used in conjunction with a transport model to 
identify the sources and estimate the emissions that produced the observed field. 
Examples of this technique include Dubovik et al., (2008), who combined MODIS AOD 
observations with the GOCART transport model to characterize global aerosol sources of 
fine particles, Yumimoto et al. (2010) who combined lidar data and model, and Wang et 
al., (2012) who combined MODIS-observed radiances with the GEOS-Chem model, 
focusing on Asian dust sources. 

Because heat from fires can propel particles quickly out of the boundary layer to higher 
layers of the atmosphere, these particles are sometimes injected into faster-flowing air 
streams. The resulting transport will be different from a model simulation that fails to 
consider such fast convective processes and instead leaves the particles in the slower 
moving boundary layer, where more efficient removal occurs. The importance of 
injection height has been discussed in context of boreal fire smoke in Canada being 
transported over the U.S. eastern seaboard (Colarco et al., 2004) and Amazon biomass 
burning smoke (Freitas et al., 2006). Satellites have made a major contribution here, 
especially the multi-angle MISR that uses stereo views and parallax to determine the 
heights of well defined smoke plumes (Kahn et al., 2007), and CALIOP that provides 
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vertical aerosol profiles (Amiridis et al., 2010). These products offer models information 
on injection heights for specific cases. The two products are complementary, as MISR 
provides greater coverage than the curtain-wide CALIOP, whereas CALIOP offers 
greater details in the vertical and provides the height of smoke after it becomes diffuse 
(Kahn et al., 2008). Val Martin et al. (2010) pushed one step beyond the individual case 
studies by establishing a climatology of injection heights using five years of MISR data 
over North America, as shown in Figure 8. This climatology li nk ed the 4-12% of the 
plumes that permeated the boundary layer, based on MISR data, with specific 
meteorological conditions and with higher fire radiative power measured by the MODIS 
sensor. Guan et al. (2010) found the linear relationship between coincident OMI AAI 
and CALIOP plume height for young biomass burning plumes, allowing the 
identification of high-altitude plumes in the multi-decadal TOMS and OMI AAI data, and 
hence making it possible to construct a data set for validating global fire plume heights in 
chemistry transport models. These associations from satellite products, begin to give the 
modeling community the information they need to develop and refine plume rises 
models, so smoke from fires can be placed at the correct altitude within the model, thus 
improving aerosol transport modeling (Sessions et al., 2011). Injection heights have been 
characterized by satellite also for dust plumes (Kahn et al., 2007; Yumimoto et al., 2009) 
and volcanic plumes (Haywood et al., 2010; Scollo et al., 2010). 

3.4.2 Model evaluation 

Satellite observations of AOD and vertical distributions have been used extensively to 
evaluate model simulations, which in many cases is the best way to bound model 
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estimates of aerosol transport. One of the earliest uses of satellites to evaluate model 
simulations of aerosol transport was Tegen et al., (1997) who showed both the 
advantages and challenges of matching model results with satellite products, in their case 
AOD from AVHRR. Their comparisons revealed both model inaccuracies and retrieval 
problems. Other studies followed, to the point that model comparison against satellite 
AOD has become almost a standard requirement. In the 1990s and 2000s, the 
comparisons relied on the TOMS and AVHRR products (e.g., Chin et al., 2002) or simple 
imagery (Colarco et al., 2003). When EOS-era aerosol satellite products became 
available, studies made use of the new, more quantitative products, mostly MODIS AOD, 
but also MISR AOD over land and ocean (Chin et al., 2004; Heald et al., 2006; Colarco 
et al., 2010; Ford and Heald, 2012; Ridley et al., 2012). Evaluation can also include 
comparison of Angstrom Exponent or size parameter (Colarco et al., 2010) and 
anthropogenic component of AOD (Yu et al., 2009b; Yu et al., 2012a). A composite of 
multiple satellite products, pre- and post-EOS era, became the standard to which the 
AEROCOM ensemble of models were held (Kinne et al., 2006). Then with the addition 
of CALIOP, the ability to evaluate the vertical distribution of aerosol became possible 
(Generoso et al., 2008; Yu et al., 2010; Koffi et al., 2012; Ridley et al., 2012). 

An evaluation of GOCART three-dimensional distributions of aerosols with MODIS and 
CALIOP measurements on a global scale is shown in Figure 9 and 10 (Yu et al., 2010). 
Figure 9 shows a comparison of GOCART AOD with collocated CALIOP and MODIS 
AOD, and Figure 10 compares the aerosol scale height between GOCART simulations 
and CALIOP measurements. The aerosol scale height is defined as a characteristic 
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altitude below which 63% of total columnar AOD is present. As shown in Figure 9, 
CALIOP, MODIS, and GOCART give generally consistent spatial patterns of AOD and 
its seasonal variations, with major continental source regions and the trans-Atlantic 
transport of Saharan dust being readily identified. Several major differences are evident 
on regional and continental scales. One of other most pronounced differences is 
associated with the intercontinental transport of aerosols. Both MODIS and GOCART 
show that the trans-Pacific transport of aerosol from East Asia to North America is fairly 
strong in MAM, with AOD greater than 0. 15 over the nearly entire mid-latitude North 
Pacific. However, the CALIOP observations show a much weaker trans-Pacific transport. 
Similar differences also exist for the trans-mid-Atlantic transport of aerosol from North 
America to West Europe. On the contrary, the westward transport of aerosol, mainly 
Saharan dust, by trade winds over tropical Atlantic is stronger and more extended from 
MODIS and CALIOP observations than the GOCART model. 

As shown in Figure 10, the GOCART aerosol scale heights are consistently higher than 
the CALIOP observations. The differences are particularly large at the polar regions and 
northern hemispheric mid-latitudes away from the source regions where aerosols are 
generally transported from outside. The long-range transport of aerosol in these regions is 
usually associated with mid-latitude cyclones that can effectively lift pollution from the 
atmospheric boundary layer (ABL) to the upper troposphere. Note that the general high 
bias of modeled aerosol scale height is not limited to GOCART model, as indicated by a 
recent AeroCom comparison that involves a dozen models (Koffi et al., 2012). As 
discussed in Yu et al. (2010), the much lower CALIOP scale height than model 
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simulations in these regions may result from the CALIOP sampling of cloud-free 
observations that may bias the scale height to low altitudes. CALIOP may miss to detect 
some optically thin layers in the FT because of the detection limit of lidar, resulting in 
lower scale heights. It is also possible that models overestimate the vertical transport of 
aerosols and gives higher scale heights. 

Comparisons of model simulations against satellite observations provide valuable hints 
for improving model representations of crucial processes such as aerosol scavenging 
removal (Bourgeois and Bey, 2011; Ridley et al., 2012) and the rise of smoke plumes 
(Sessions et al., 2011; Dirksen et al., 2009). Most recently, Ridley et al., (2012) 
demonstrate the extent to which a suite of satellite products can be used to evaluate a 
model’s ability to simulate aerosol transport, in this case multi-year Saharan dust 
transport over the Atlantic. The model is shown to capture the variability of these dust 
plumes by correlation with MODIS and MISR and to match the vertical distribution of 
the plumes as they exit the African continent by quantitative comparison with CALIOP. 
They also find that the lifetime of the simulated dust is a few days shorter than that 
observed from satellite. Comparison of model rain rates with TRMM satellite-derived 
rain rates shows that wet deposition in the model is too strong. 

3.4.3 Integrated approach 

Integrating satellite measurements and aerosol transport models is necessary for a 
comprehensive characterization of aerosol intercontinental transport. On one hand, 
satellite measurements provide observational evidence of aerosol intercontinental 
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transport that can be quantitative (Section 3.3). However, these observations are 
restricted by limited satellite capabilities such as observing only once or twice daily, and 
for most sensors, only under cloud-free conditions. Space-bome lidars add an essential 
third dimension to the characterization, but lidars are limited by extremely narrow swaths 
and poor spatial sampling. On the other hand, model simulations offer a much more 
comprehensive picture of aerosol intercontinental transport by providing multi- 
dimensional characteristics of aerosol in a continuous manner, including clearly 
distinguishing by aerosol species. But they are limited by their imperfect ability to 
simulate physical and chemical processes in the atmosphere. A number of studies have 
used this satellite-model integrated analysis approach to characterize the intercontinental 
transport of dust and pollution aerosol (Eguchi et al., 2009; Hara et al., 2009; Itahashi et 
al., 2010; Uno et al., 2008, 2009, 2011; Yumimoto et al., 2009, 2010). 

Satellite aerosol products can be assimilated into run-time models to actively nudge 
models towards real-world conditions, often producing better representation of the 
aerosol field and aerosol transport. This is especially useful for real-time aerosol 
forecasting. One of the earliest examples of aerosol assimilation involved the Model of 
Atmospheric Transport and Chemistry (MATCH) and how it was applied to the INDOEX 
experiment in 1999 (Collins et al., 2001). The model was run to produce an aerosol field, 
then that field was adjusted to better resemble the AOD as derived from AVHRR satellite 
observations. The results, both model-alone and model-with-satellite assimilation, were 
compared to shipboard sunphotometers; the assimilated results show much closer 
agreement to the sunpho tome ter measurements. The assimilation provided aerosol 
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information when satellite retrievals were impossible: at night and in cloudy situations 
(Rasch et ah, 2001). In recent years, other modeling groups have assimilated EOS-era 
satellite aerosol information for real-time aerosol forecasting. MODIS aerosol optical 
depth, over land and ocean, is automatically ingested and assimilated into the European 
Centre for Medium Range Forecasting (ECMWF) for aerosol predictions (Benedetti et 
al., 2009). Figure 11 shows comparisons between AODs from the free-running ECMWF 
model and the ECMWF-MODIS analysis compared to MODIS observations for 5-6 
March 2004 Saharan dust outbreak. Figure 12 shows a corresponding comparison 
between AOD at 670 mn from the analysis and the free-running forecast and AERONET 
AODs at 675 mn for Agoufou (Mali), Dakar (Senegal) and Cape Verde. Although the 
satellite observed shape of the dust outflow is well represented in the free-running model, 
the magnitude of the AODs is much lower than both MODIS and AERONET 
measurements. The analysis shows much better agreement with AERONET 
measurements. The peaks shown in the AERONET data were well captured by the 
analysis, with the exception of the 8-9 March AOD maximum over Agoufou. 

MODIS aerosol optical depth is also assimilated into the Naval Research Laboratory 
(NRL) Atmospheric Variational Data Assimilation System (NAVDAS) over ocean 
(Zhang et al., 2008) and is being prepared for assimilation over land (Hyer et al., 2011). 
Likewise, CALIOP vertical profiles (Campbell et al., 2010) and MISR-derived aerosol 
properties have also been assimilated into the same system. By adding these satellite data 
products, model forecasts and aerosol transport results were quantitatively improved 
(Zhang et al., 2011). Other assimilation studies include Yumimoto and Takemura (2011) 
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who assimilated one month of Terra and Aqua MODIS AOD. Studies have also used 
satellite products to constrain model simulations through diagnostic analysis (Yu et al., 
2003; Matsui et al., 2004; Liu et al., 2005). However data assimilation and analysis 
currently provide more total AOD than component AOD constraints, due to the lack of 
accurate, global aerosol type observations from satellites. 

Satellite measurements can also be integrated with model systems to assess the influences 
of aerosol intercontinental transport on atmospheric composition, air quality, weather, 
and climate. Over ocean, the intercontinental transport of continental aerosols can often 
be distinguished from relatively clean marine aerosols and the interactions of aerosols 
with cloud and precipitation may be inferred from a correlative analysis of satellite 
measured aerosols, clouds, and precipitation (e.g., Kaufman et al., 2005d; Koren et al., 
2005; Zhang et al., 2007), though significant uncertainty still apply. Over land, however, 
distinguishing ICT aerosols from local sources remains a great challenge for satellite 
remote sensing. Aerosols from different continents usually have similar characteristics; 
and current satellites are not capable of discerning the usually marginal differences 
between continents. Thus assessing the impacts of aerosol ICT on air quality, weather, 
and climate in downwind continents has to rely on model systems that link important 
atmospheric processes at local, regional, and hemispheric scales. Global models are 
needed to provide regional models with appropriate boundary and initial conditions, 
whereas regional models with higher resolution are needed to adequately represent major 
chemical and meteorological processes important for air quality, weather, and climate. 
Considering large uncertainties associated with global aerosol models, measurements 
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from satellites should be used to provide chemical boundary conditions of aerosol ICT 
for regional models. 

4. Conclusions and Outlooks 

Satellites provide ideal platfonns for studying aerosol intercontinental transport, because 
of the frequent sampling over global and multi-year scales. In the last decade or two, 
major advances have been achieved in both passive and active aerosol remote sensing 
from space, providing better coverage, higher spatial resolution, improved AOD 
accuracy, emerging climatology of vertical distributions, and enhanced capabilities to 
characterize aerosol type, as provided by observational constraints on particle properties. 
Satellites also provide correlative measurements of trace gases, as well as the strengths 
and properties of aerosol source. These advances offer unprecedented opportunities to 
develop a measurement-based assessment of aerosol intercontinental transport. 

The accumulation of advanced satellite observations over the past decade has 
significantly improved the understanding of aerosol intercontinental transport and its 
impacts on climate and air quality, from largely qualitative characterization to 
increasingly quantitative assessment. The observations have been used to characterize the 
three-dimensional evolution of aerosol loading and particle properties (size, and shape) 
during intercontinental transport on an episodic basis. Long-term measurements with 
daily frequency have been used to characterize seasonal and inter-annual variations of 
aerosol ICT and associated controlling factors. Most recently, estimates of aerosol export 
and import have been performed using the satellite measurements of aerosol optical 
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depth, vertical distribution, and microphysical properties over ocean. The findings have 
important implications for air pollution control and climate mitigation. 

Measurement-based quantitative assessment of aerosol ICT requires a synthesis of 
multiple data sets. No sensor, known or expected in the foreseeable future, has sufficient 
capability to fully characterize aerosol ICT without additional infonnation. For example, 
MODIS offers among the most complete coverage, which is essential to capturing 
specific transport events, but cannot directly identify transport height or aerosol type. 
MISR provides much better characterization of aerosol particle properties and how those 
evolve along the transport path, but has a more limited sampling. CALIOP is the ideal 
instrument to identify transport height and can distinguish dust from combustion aerosol, 
but is extremely limited in its sampling. The other satellite sensors also each make their 
own specific contributions. Combining information from different sensors is shown to 
provide quantitative estimates of the material flux imported and exported at continental 
boundaries, albeit still with large error bars. Narrowing these error bars requires 
continued characterization of the aerosol using additional data sets and new validated 
products from the current suite of satellite sensors (i.e., more accurate particle properties, 
AOD above clouds). In turn, validation of the satellite data sets requires a concerted 
effort from suborbital observations. 

In addition to the satellite-alone observation-based applications, satellite measurements 
have also provided several critical inputs for constraining chemical transport models, and 
valuable datasets for evaluating and improving model simulations of aerosol 
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intercontinental transport. The ultimate means to quantify aerosol ICT will necessarily 
involve a combination of multiple satellite observations and models. This may take the 
form of directly assimilating satellite aerosol products into models, or constraining 
models with observational results. These constraints amount to more than just requiring 
the model to match the satellite-derived AOD, or even the satellite-derived particle 
properties. The models also need to simulate the satellite-derived material fluxes into and 
out of continents at the continental boundaries. It is recommended that AeroCom and 
HTAP model experiments consider including the export/import flux requirement into 
their data protocols, so that models can be evaluated with satellite observations. 

Despite this significant progress, substantial effort is needed in the future to fully explore 
the potential offered by existing satellites, and to develop new and enhanced satellite 
missions: 

(a) Geostationary satellite measurements have been underused. Although such 
measurements usually do not have adequate spatial coverage to track the 
evolution of aerosol plumes along the entire intercontinental transport path, in 
combination with polar orbiting satellite measurements they can provide the hour- 
to-hour evolution of the plumes in many regions of interest. 

(b) Satellite-based studies should extend to all major ICT paths. The trans-Atlantic 
transport of Saharan dust has been well studied, and there are studies addressing 
the trans-Pacific transport of dust and pollution. However, except for a few 
efforts, the Indian Ocean, the Southern Hemisphere, pollution transport from 
North America has not been addressed by observationally based work. 
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(c) More information can be extracted from existing sensors. The integrated use of 
multiple satellite measurements can potentially improve the characterization of 
aerosol temporal and spatial distribution, distinguish aerosol types and overcome 
uncertainties introduced by sampling issues. In particular, integration of non- 
aerosol products such as carbon monoxide and SO 2 gas can add valuable 
information, but are mostly missing from current studies. Novel analysis 
techniques that include combining absorption Angstrom exponent and extinction 
Angstrom exponent for better aerosol type characterization (Russell et al., 2010) 
or retrieving AOD from above clouds using passive sensors (Waquet et al., 2009; 
Torres et al., 2012; Yu et al., 2012b) could also aid in improving assessment of 
aerosol ICT from satellite measurements alone. 

(d) Satellite products require suborbital validation. As more infonnation is extracted 
from existing sensors, these new products require validation that only suborbital 
measurements can supply. Validation has to go beyond comparing AOD with 
AERONET, to make use of the detailed particle property measurements being 
made from airborne and ground-based instruments. 

(e) Continuity and enhancement of satellite capabilities. Now that we are fully 
exploiting the information content and synergy of the EOS-era satellites, we note 
that these sensors are aging and are operating almost a decade beyond their initial 
mission lifetimes. New and future satellites (ie. Suomi-NPP VIIRS, GOES-R 
ABI, the Earthcare suite of sensors) can continue some, but far from all EOS-era 
capability. Key to future abilities is the Earthcare ATLID, a lidar to succeed 
CALIOP, although its different spectral configuration raises questions of 
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continuity. We are missing a next-generation multi-angle, polarimeter imager that 
can provide the detailed particle properties we have come to expect from MISR 
and POLDER. 

(f) Consequences of aerosol ICT should be explored. What is the radiative effect of 
this transport over oceans and over continents? What is the effect on air quality? 
How are these consequences partitioned between dust and combustion aerosols? 
How much is anthropogenic and how much is natural? The answers to these 
questions will require a synthesized combination of multiple observation systems 
and multi-scale models. However, satellite -based observations can take a leading 
role, and make significant quantitative contribution to finding the answer to these 
important questions. 

(g) Satellite-based assessments alone are insufficient to answer all the questions. 
Modeling studies will be necessary, but models must use the information gained 
by satellite studies to constrain emissions, evaluate AOD distributions and 
compare with import/export assessments at continental boundaries. Assimilation 
provides one way forward, but sacrifices the ability of the model studies to 
understand aerosol processes in the atmosphere through the exercise of 
modeling/evaluation/model improvement. 
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Acronyms 

AAI: 

Absorbing Aerosol Index (sometimes as AI) 

A AT SR: 

Advanced Along-Track Scanning Radiometer 

ABI: 

Advanced Baseline Imager 

ABL: 

Atmospheric Boundary Layer 

ACA: 

Above Cloud Aerosol 

ACE: 

Aerosol-Cloud-Ecosystem (NASA mission in planning) 

ADEOS: 

Advanced Earth Observing Satellite 

AE: 

Angstrom Exponent 

AEROCOM: 

Aerosol Comparisons between Observations and Models 

AERONET: 

AErosol RObotic NETwork 

AIRS: 

Atmospheric Infrared Sounder 

AOD: 

Aerosol Optical Depth 

ATLID: 

Atmospheric Lidar 

AT SR: 

Along Track Scanning Radiometer 

AVHRR: 

Advanced Very High Resolution Radiometer 

BUV: 

Backscatter UV 

CALIOP: 

Cloud- Aerosol Lidar with Orthogonal Polarisation 

CALIPSO: 

Cloud-Aerosol Lidar with Orthogonal Polarisation 

CLAES: 

Cryogenic Limb Array Etalon Spectrometer 

CTM: 

Chemical Transport Model 

ECMWF : 

European Centre for Medium Range Forecasting 

EOS: 

Earth Observing System 

EP: 

Earth Probe 

ERS: 

European Remote Sensing 

FMF: 

Fine Mode Fraction (aerosol) 

FT: 

Free Troposphere 

GEOS-Chem: 

Goddard Earth Observing System (GEOS) - Chemistry model 

GLAS: 

Geoscience Laser Altimeter System 

GMI: 

Global Modeling Initiative 

GOCART : 

Goddard Chemistry Aerosol Radiation and Transport model 

GOES: 

Geostationary Operational Environmental Satellites 

GOME: 

Global Ozone Monitoring Experiment 

HSRL: 

High Spectral Resolution Lidar. 

HTAP: 

Hemispheric Transport of Air Pollutants 

HYSPLIT: 

Hybrid Single Particle Lagrangian Integrated Trajectory Model 

IASI: 

Infrared Atmospheric Sounding Interferometer 

ICT: 

InterContinental Transport 

INDOEX: 

Indian Ocean Experiment 

ITCZ: 

Inter-Tropical Convergence Zone 

LH: 

Layer Height (aerosol) 

LITE; 

Lidar In-Space Technology Experiment 
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MAI AC: 

MultiAngle Implementation of Atmospheric Correction 

MAN: 

Marine Aerosol Network 

MATCH: 

Model of Atmospheric Transport and Chemistry 

MERIS: 

MEdium Resolution Imaging Spectrometer 

MISR: 

Multi-angle Imaging SpectroRadiometer 

MJO: 

Madden- Julian Oscillation 

MLS: 

Microwave Limb Sounder 

MODIS: 

Moderate Resolution Imaging Spectroradiometer 

MOPIIT : 

Measurements Of Pollution In The Troposphere 

MSG: 

Meteosat Second Generation 

NAO: 

North Atlantic Oscillation 

NAVDAS: 

NRL Atmospheric Variational Data Assimilation System 

NPP: 

NPOESS Preparatory Project 

NRC: 

National Research Council 

NRL: 

Naval Research Laboratory 

NSF: 

Non-Spherical Fraction (aerosol) 

OMI: 

Ozone Monitoring Instrument 

PARASOL: 

Polarization and Anisotropy of Reflectances for Atmospheric Sciences 
Coupled with Observations from a Lidar 

PM: 

Particulate Matter 

POLDER: 

Polarization and Directionality of the Earth’s Reflectance 

RAIR: 

Relative Annual Intercontinental Response 

RH: 

Relative Humidity 

SAGE: 

Stratospheric Aerosol and Gas Experiment 

SAM: 

Sample Analysis at Mars 

SCIAMACHY : 

Scanning Imaging Absorption Spectrometer for Atmospheric 
Chartography/Chemistry 

SeaWiFS: 

Sea-viewing Wide Field-of-view Sensor 

SEVIRI: 

Spinning Enhanced Visible Infra-Red Imager 

SSA: 

Single-Scattering Albedo (aerosol) 

TES: 

Thermal Emission Spectrometer 

TOA: 

Top-Of- Atmosphere 

TOMS: 

Total Ozone Mapping Spectrometer 

TRMM: 

Tropical Rainfall Measuring Mission 

UV: 

Ultra-violet (radiation) 

VIIRS: 

Visible Infrared Imaging Radiometer Suite 

VP: 

Vertical Profile (aerosol) 

WCB: 

Warm Convey Belt 
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Table 1: Summary of major operational and research products from satellite remote sensing that can be applied 
for characterizing aerosol intercontinental transport. 


Sensor/platform 

Swath 

Measurement 

Period 

Column 

load 

Particle Properties 
(for aerosol type characterization) 

Transport Heights 

AOD 

AE 

FMF 

/size 

SSA 

AAI 

NSF 

LH 

VP 

ACA 

TOMS/Nimbus, 
ADEOS, EP 

-2700 km 

1979 -2001 

V 




V 




V 

AVHRR series 

-2400 km 

1981 - present 

V 

(O) 

V 

(O) 








A(A)TSR/Envisat 

-500 km 

1995 - present 

V 

V 








GOME/ERS-2 

GOME-2/METOP 

-2000 km 

1995 -2003 
2007 - present 

V 




V 





POLDER/ADEOS 
1, 2/PARASOL 

-2400 km 

1996, 2003 
2004 - present 

V 

V 

V 

V 


V 



V 

SeaWiFS/SeaStar 

-2800 km 

1997 -2010 

V 

V 








MISR/Terra 

-360 km 

2000 - present 

V 

V 

V 

V 


V 

V 



MODIS/Terra 

MODIS/Aqua 

-2300 km 

2000 - present 
2002 - present 

V 

V 

V 

(O) 

V 

(L) 






MERIS/Envisat 

-1150 km 

2003 - present 

a/ 

a/ 








SCIAMACHY/ 

Envisat 

-960 km 

2003 - present 

V 




V 





AIRS/Aqua 

-1600 km 
(night) 

2004 - present 

a/ 

(dust) 






a/ 

(dust 

) 



OMI/Aura 

-2600 km 

2004 - present 

V 



V 

V 




V 

GLAS/ICESat 

-70 m, day & 
night 

2003 -present 

V 







V 


CALIOP/CALIPS 

0 

-70 m, day & 
night 

2006 - present 

V 

V 




V 



V 

IASI/MetOp 

-2050 km 

2006 - present 

a/ 

(dust) 






V 

(dust 

) 



VIIRS/NPP 

Soumi 

-3000 km 

2012 - present 

V 

V 

(O) 

V 

(O) 







VISSR/GOES 

series 

-half-hourly, 
North and 
Central 
America 

2003 - present 

V 









SEVIRI/MSG 

-hourly, 

West Europe, 
Africa, and 
the Atlantic 

2004 - present 

V 










Acronyms: L - Land only; O - Ocean only; AOD - aerosol optical depth; AE - Angstrom exponent; FMF - fine- 
mode fraction (with respect to AOD); SSA - single-scattering albedo; AAI - absorbing aerosol index; NSF - non- 
spherical fraction of AOD; LH - layer height; VP - vertical profile; ACA - above-cloud aerosol. 
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Table 2: Comparisons of trans-Pacilic dust transport efficiency (a ratio of dust mass flux in North 
America inflow to that in East Asia outflow) and dust deposition to the North Pacific. 


Transport 

efficiency 

Dust 

deposition 
(Tg a 1 ) 

Region 

Notes 

Source 

Satellite-based estimate 

40% 

84 

30°-60°N, 

140°E-130°W 

Deposition estimate should bias 
high because of the existence of 
northward transport to Arctic 

Yu et al. (2012c) 

Estimate based on in-situ measurements 

- 

10 

East China 
Sea only 

Using 20g nf 2 a' 1 deposition rate 
over the 0.5xl0 12 m 2 area 

Hsu et al. (2009) 

- 

67 

(18 -260) 

East China 
Sea only 

Estimated based on surface 
concentration measurements and 
model-based deposition rates 

Gao et al. (1997) 

- 

96 

North Pacific 

Prospero et al. (1996) 

- 

480 

North Pacific 

Duce et al. (1991) 

- 

20 

0°-50°N, 

150°E-130°W 

Uematsu et al. (1985) 

Model simulations 

34% 

88 

30°-60°N, 

140°E-130°W 

GOCART model; including 
northward transport to Arctic 

Yu et al. (2012c) 

25% 

92 

30°-60°N, 

140°E-130°W 

GMI model; including northward 
transport to Arctic 

Yu et al. (2012c) 

40% 

- 

25°-55°N, 

150°E-130°W 

2007 

Su and Toon (2011) 

15% 

- 

120°E-120°W 

For events in mid- August of 2009 

Yumimoto et al. 
(2010) 

30% 

- 

0°-70°N, 

140°E-120°W 

For events during May 5-15, 2007 

Eguchi et al. (2009) 

43% 

- 

0°-90°N, 

120°E-120°W 

Highly lofted (8- 10km) dust event 
in May 2007 

Uno et al. (2009) 

35% 

- 

140°E-130°W 

Taklimakan dust storm in May 
2007 

Yumimoto et al. 
(2009) 

- 

68 

North Pacific 

3-model composite; constrained 
by observations 

Jickells et al. (2005) 

- 

56 

North Pacific 


Tegen et al. (2004) 

- 

35 

North Pacific 


Luo et al. (2003) 

- 

64 

western 
North Pacific 
(-25% of the 
basin) 

March 1994-May 1995 

Uematsu et al. (2003) 

- 

31 

North Pacific 


Zender et al. (2003) 

- 

92 

North Pacific 


Ginoux et al. (2001) 
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Figure 1: Multi-model derived average relative annual intercontinental response (RAIR) 
for surface concentration of S04+POM and aerosol optical depth (AOD) by chemical 
component in four major pollution regions (NA - North America, EU - Europe, EA - 
East Asia, SA - South Asia), which is based on HTAP source/receptor relationship 
experiments by reducing anthropogenic emissions by 20% in each of the four regions. 
RAIR in a receptor region represents the percentage contribution of the intercontinental 
transport of foreign sources relative to the sum of foreign and domestic sources. Standard 
deviations, as indicated by error bars, reflect the variability of the 10 models in the 
ensemble. Figures are adapted from HTAP (2010) for surface concentration and Yu et 
al. (2012a) for AOD. 
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Figure 2: Spatial variations of aerosol type on April 11, 2001 as revealed by a composite of aerosol 
optical depth (at 550 nm) and fine-mode fraction observed by MODIS/Terra (top panel) and simulated 
by GOCART model (bottom panel) for April 11, 2001. Industrial pollution and biomass burning 
aerosols are predominated by small particles (shown as red), while mineral dust consists of a large 
fraction of large particles (shown as green). Bright red and bright green indicate heavy pollution and 
dust plumes, respectively (image credit: Reto Stockli and Yoram Kaufman). 
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Figure 3. MISR oblique-camera images and retrievals of optical depth, Angstrom exponent, and non- 
spherical fraction of aerosol plumes. Top: southern Mexico smoke plumes drifting northward over the 
Gulf of Mexico on May 2, 2002 (Terra orbit 12616) Bottom: dust blowing off Sahara desert to Cape 
Verde islands on February 6, 2004 (Terra orbit 22006). The width of each image it the full MISR 
swath - about 400 km acorss; the vertical dimension is presented at the same spatial scale. (. Adapted 
from Diner et al., 2005. Image credit: MISR Team, Jet Propulsion Lab/Caltech and NASA Goddard 
Space Flight Center) 
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Figure 4a: An A-Train integrated characterization of dust and smoke aerosols transported over the 
Gulf of Guinea on January 31, 2008. Red line in maps of aerosol and CO horizontal distribution 
indicates the CALIPSO track. 
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CDC1 Meteorological Data 



Figure 4b: 6-day HYSPLIT backward trajectories ending at [9°N, 15°W] in the smoke-dominated 
regime and at [4°N, 7.5°W] in the dust-smoke mixed regime on January 31, 2008 1400UTC. The air 
mass arrived at both locations passed over the smoke-active Sahel region during the January 28-30 
period. During the January 26-27 period, the air mass arriving at [4°N, 7.5°W] originated from the 
Bodele Depression - the major source of Saharan dust in the winter, whereas that arriving at [9°N, 
15°W] drifted off north to the Bodele Depression and were influenced less by dust. 
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Figure 5. A three-dimensional view of the hemispherical transport of Asian dust during May 9-22, 
2007 as provided by CALIOP measurements of vertical distributions of depolarization ratio (a, dust is 
distinguished from other types of aerosols by exhibiting high depolarization ratio) and extinction 
coefficient (b), and a global model simulation of extinction coefficient in the dust layer (c). The 
vertical and horizontal transport paths are further illustrated by HYSPLIT trajectories. (Figure taken 
from Uno et al., 2009, with permission from Nature Geoscience ) 
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Figure 6: Seasonal variations of trans-Atlantic transport of dust AOD observed by MODIS, MISR, 
and PARASOL, and simulated by GOCART model. 
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Figure 7: Meridional variations of annual trans- 
pacific aerosol import to North America as 
estimated from MODIS and CALIOP 
measurements (a) and simulated by GOCART 
(b) and GMI (c) model, which is summarized 
from results documented in Yu et al. (2008) and 
Yu et al. (2012c). 
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Figure 8. MISR observed smoke plume heights (above the ground level) in each biome over North 
America in 2002, and 2004-2007. Bar plots indicate the distribution of the data. The medians (red 
circles) and the means (black squares) are shown along with the central 67% (color coded box) and the 
central 95% (thin black lines). Data that fall outside the bar plots are plotted with black circles. (Figure 
taken from Val Martin et al., 2010, with permission from Atmospheric Chemistry and Physics). 
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Figure 9. Comparisons of seasonal AOD between GOCART model, CALIOP and MODIS retrievals. 
( adapted from Yu et al., 2010, with permission from American Geophysical Union) 
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Figure 10: Comparison of aerosol scale height between GOCART model and CALIOP measurements 
( adapted from Yu et al., 2010, with permission from American Geophysical Union). 
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Figure 11. AOD (at 550 nm) distributions of March 2004 Saharan dust outbreak as simulated by the 
ECMWF model: (a) free-running model, (b) assimilation with MODIS observations, and (c) MODIS 
observations for (left) 5 March 2004 at 1200 UTC and (right) 6 March 2004 at 1200 UTC. (Figure 
taken from Ben edetti et al., 2009, with permission from American Geophysical Union). 
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c Comparison of exlz & ezub AOT at 670 nm and L2.0 Aeronet AOT at 675 nm over 

Capo_Verde (lat = 1 6.73, Ion = -22.93). Period = 28/02/2004 - 1 5/03/2004. FC start hrs = 0, 1 2Z 



Figure 12. Comparisons of ECMWF analysis AOD at 675 nm with AERONET observations for the 
Saharan dust outbreak of March 2004 at: (a) Agoufou (Mali), (b) Dakar (Senegal), and (c) Cape Verde. 
AERONET data are shown in light blue, the analysis is in red, and the free-running forecast is in dark 
yellow. (Figure taken from Benedetti et al., 2009, with permission from American Geophysical Union). 
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Supplementary Online Materials 
for 

Satellite Perspective of Aerosol Intercontinental Transport 
From Qualitative Tracking to Quantitative Characterization 

Hongbin Yu, Lorraine A. Remer, Ralph A. Kahn, Mian Chin, Yan Zhang 

A. Descriptions of Major EOS-era Aerosol Sensors 

NASA’s Earth Observing System (EOS) is an ambitious program designed in the 1980s, implemented 
in the 1990s and executed beginning with the launch of the Terra satellite at the end of 1999. The 
program has resulted in a constellation of polar orbiting satellites, equipped with sensors that provide 
observations in support of Earth system science. The main satellites and sensors providing aerosol 
information include: Terra-MODIS, Terra-MISR, Aqua-MODIS, Aura-OMI, PARASOL-POLDER 
and CALIPSO-CALIOP. 

A.l. MODIS multi-wavelength measurements: Each MODIS sensor obtains near global, daily 
observations of atmospheric aerosols with a swath width of 2330 km and 36 channels ranging from 
0.41 to 15 pm. Eight channels, nominally, 0.41, 0.47, 0.55, 0.66, 0.87, 1.24, 1.63 and 2.13 pm are used 
to retrieve aerosol properties in cloud-free scenes with appropriate surface features (Martins et al., 
2002; Li et al., 2005; Remer et al., 2005; Hsu et al., 2006). Because of its wide spectral range, and the 
simplicity of the dark ocean surface, MODIS has the capability of retrieving aerosol optical depth with 
a relatively high accuracy of ±0.03±0.05x (Remer et al., 2002; 2005; Levy et al., 2010). The ocean 
conditions also pennit retrieval of information on particle size. This can take the fonn of Angstrom 
Exponent, effective radius or fine mode fraction (See Eq. 3 of Section 3.2). The three retrieved size 
parameters are all inter-related and represent the same basic information content available from the 
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spectral signature of the satellite-measured radiances against the dark ocean surfaces. Over vegetated 
land, the MODIS dark-target algorithm estimates surface reflectance in the blue and red channels 
based on an empirical relationship using the shortwave infrared radiance and retrieves AOD (Kaufman 
et al., 1997; Levy et al., 2007). AOD is retrieved at 0.55 pm with an accuracy of ±0.05±0.15x and then 
extrapolated to 0.47 pm and 0.66 pm using the selected aerosol model (Chu et al., 2002; Remer et al., 
2005; Levy et al., 2010). Because of the complex spectral structure of the land surface (note that 
vegetation is bright, not dark in the near-IR), retrieval of aerosol information beyond AOD is not 
possible. 

Alternative aerosol algorithms have been developed and applied to MODIS radiances. The Deep Blue 
(Hsu et al., 2004) and the MultiAngle Implementation of Atmospheric Correction (MAIAC) 
algorithms (Lyapustin et al., 2011) offer complementary and alternative products to the original Dark 
Target algorithm described above. The Deep Blue algorithm takes advantage of the much darker 
surface reflectance in the blue than at red wavelengths (Hsu et al., 2004). This is especially useful for 
retrieving AOD over bright surfaces such as deserts, although Deep Blue is not limited to desert 
scenes. MAIAC retrieves aerosol information over land simultaneously with parameters of a surface 
bidirectional reflectance factor model using image-based processing applied to the time series of 
MODIS measurements (Lyapustin et al., 2011). The Deep Blue products are being produced 
operationally from Aqua-MODIS data (currently not Terra-MODIS). The MAIAC products are not 
operational at this time. There is a possibility that the time series of MODIS aerosol measurements 
may be continued from measurements of the Visible/Infrared Imager Radiometer Suite (VIIRS) 
onboard the most recently launched Suomi National Polar-orbiting Partnership (NPP). VIIRS has 
aerosol remote sensing capabilities similar to those of MODIS (Lee et al., 2006), with a broader swath 
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width than MODIS that fully covers the Earth every day. However, the spectral response of each 
channel is not identical to MODIS and other differences require analysis to determine to what accuracy 
can VIIRS continue the MODIS aerosol record. 

A.2. MISR multi-angle measurements: MISR, aboard the sun-synchronous polar orbiting Terra 
satellite, measures upwelling solar radiance in four spectral bands (centered at 446, 558, 672, and 866 
mn) at each of nine view angles spread out in the forward and aft directions along the flight path (at 
nadir, ±70.5°, ±60.0°, ±45.6°, and ±26.1° of nadir) (Diner et al., 2002). It acquires global coverage 
about once per week. Having observations taken through air-mass-factors ranging systematically from 
one to three gives the instrument greater sensitivity to optically thin aerosol. MISR’s wide range of 
along-track view angles also makes it feasible to more accurately evaluate the surface contribution to 
the TOA radiances and hence retrieve aerosols over both ocean and land surfaces, including bright 
desert aerosol source regions and ocean regions that are contaminated by sun glint for mono- 
directional instruments (Diner et al., 1998; Martonchik et al., 1998; 2002; 2009; Kahn et al., 2005; 
2010; Kalashnikova and Kahn, 2006; 2008). MISR constraints on aerosol type also enhance MISR 
AOD accuracy compared to single-angle instruments in many circumstances. MISR provides particle 
type classification, based on particle size, shape and SSA constraints. These amount to distinguishing 
three-to-five particle size bins, two-to-four bins in SSA, and spherical vs. non-spherical particles, 
provided the mid-visible AOD exceeds 0.15 or 0.2, and the scene meets other basic retrieval-quality 
criteria (Kahn et al., 2001; 2010). Cloud screening includes tests based on location-and-time specific 
camera-by-camera radiometric thresholds, reflecting-layer stereo-derived height above the surface; and 
angular smoothness and correlation (Martonchik et al., 2009). Stereo-derived reflecting-layer heights 
also provide maps of aerosol plume elevation near sources, where contrast features can be matched in 
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multi-angular views, at 1.1 km horizontal resolution and about 0.5 km vertical resolution. More details 
about the MISR aerosol products are given in the references cited above. 

A.3. POLDER multi-directional and polarization measurements: There have been three POLDER 
instruments launched into space ADEOS-POLDER (1996-1997), ADEOS2-POLDER (2003) and 
PARASOL-POLDER (2005 - still operating at time of writing). PARASOL-POLDER flew in 
formation with the other A-Train sensors until 2010, when fuel considerations caused the satellite to 
drop to a lower orbit. Each POLDER sensor consists of a wide field-of-view imaging spectro- 
radiometer capable of measuring multi-spectral, multi-directional, and polarized radiances with a 
global coverage within two days (Tanre et al., 2011). POLDER measures in nine spectral channels 
from 0.443 pm to 1.020 pm, and in several polarization orientations at 0.49 pm, 0.67 pm and 0.87 pm. 
As the satellite passes over a surface target that target will be viewed at different spots in the image, at 
different angles. This provides multi-angle observations for each target, ±51° along track. 

The observed multi-angle polarized radiances can be exploited to better separate atmospheric from 
surface contributions over both land and ocean (Deuze et al., 2001). Over ocean, the total and polarized 
radiances at 670 and 865 mn are used to retrieve total AOD (Deuze et al., 1999), with a typical 
accuracy of ±0.05±0.05x (Goloub et al., 1999). When the geometrical conditions are optimal 
(scattering angle ranging between 90°-160°), POLDER can derive the shape of coarse-mode particles 
(Herman et al., 2005). Over land, the aerosol retrieval is based on measurements using only polarized 
light at the two wavelengths (Hennan M. et al., 1997; Deuze et al., 2001), to capitalize on the small 
and fairly spectrally neutral polarized reflectance typical of land surfaces. Because large aerosol 
particles polarize sunlight much less than sub-micron particles (except in very intense dust storm 
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events), the POLDER over-land AOD retrieval is largely sensitive to spherical, accumulation-mode 
particles (Tanre et al., 2011), providing the fine mode aerosol optical depth but not the total. 

Recently, advanced algorithms are being developed that will further exploit the infonnation content of 
the POLDER measurements (Dubovik et al., 2011; Kokhanovsky et al., 2010). The retrieval is 
designed as a statistically optimized multi-variable fitting of the complete PARASOL observation set: 
all wavelengths, all angles, and all states of polarization. The algorithm allows for a large number of 
unknown parameters and retrieves a set of parameters affecting measured radiation. The algorithm 
simultaneously retrieves aerosol and surface properties and should provide infonnation about aerosol 
sizes, shape, absorption and refractive index. 

A.4. OMI near-UV measurements: OMI flies on the Aura satellite with a swath width of 2600 km that 
covers the Earth in one day. OMI is a hyperspectral instrument that measures from 0.27 to 0.50 pm in 
bandwidths of 1 pm or less. The Aerosol Index (Al; sometimes refened to as the UV Aerosol Index, 
UVAI; or the Absorbing Aerosol Index, AAI) provides a measure of absorbing aerosol and a simple 
index to distinguish between absorbing and non-absorbing aerosols elevated above the ABL. AAI is 
calculated without assuming aerosol particle properties. It has been frequently used to characterize 
aerosol transport. 

There are two operational aerosol retrievals making use of OMI measurements, the OMAERUV and 
OMAERO. The OMAERUV algorithm uses the measured radiances at 0.354 and 0.388 pm to derive 
AOD and absorbing AOD at 0.388 pm (Torres et al., 2007). Because the surface is obscured by 
atmospheric gas Rayleigh scattering at these wavelengths, and in addition, because most surfaces have 
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very low and homogeneous reflectance in the UV range, OMAERUV retrievals cover global ocean and 
land free of snow and ice. The algorithm is based on the deviation of observed radiances from 
expected values from Rayleigh scattering and is sensitive to AOD, aerosol absorption and height of the 
aerosol layer. With only two pieces of information from the two measured channels, OMAERUV 
cannot retrieve all three unknowns. The operational algorithm assumes aerosol layer height from 
climatology, specifically the climatology as observed from CALIOP lidar, and returns AOD in the UV 
and infonnation on aerosol absorption. The AOD is extrapolated to the visible regions using the 
selected aerosol spectral model (Torres et al., 2002, 2007). The OMAERO algorithm uses up to 19 
channels in the spectral range of 0.33 pm to 0.50 pm to derive quantitative aerosol information. The 
oxygen absorption band at 0.477 pm is used to enhance the sensitivity to the aerosol layer height 
(Veihehnann et al., 2007) in a way similar to its use in the retrieval of cloud height (Acarreta et al., 
2004), and becomes one of the retrieved parameters. The primary retrieved parameters are AOD and 
the best fitting aerosol model. Both OMAERUV and OMAERO retrieve quantitative AOD. This is 
AOD of the entire column, including the surface layer. 

Due to its low spatial resolution, i.e., 13 x 24 km 2 at nadir, OMI AOD retrievals are prone to sub-pixel 
cloud contamination. OMI AOD values are generally biased high with respect to MODIS 
measurements, likely as a result of a calibration offset and subpixel cloud contamination (Ahn et al., 
2008). In comparison, OMI absorbing AOD is less sensitive to the cloud contamination because of a 
partial cancellation of cloud contamination effects on AOD and single-scattering albedo (Torres et al., 
2007). 
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A.5. CALIOP vertical profiling: Following the Lidar In-space Technology Experiment (LITE) 
demonstration mission in 1994 (Winker et ah, 1996), the Geoscience Laser Altimeter System (GLAS) 
was launched in early 2003, becoming the first polar-orbiting satellite lidar (Schulz, 2001; Zwally et 
ah, 2002). It detects clouds and aerosols (Spinhirne et al., 2005) for one month out of every three-to- 
six, due to the unfortunate failure of a laser. Because of its inconsistency of making measurements, 
GLAS has not often been used for aerosol transport studies. 

In 2006, the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) was 
launched, carrying a two-wavelength, polarization lidar, namely the Cloud-Aerosol Lidar with 
Orthogonal Polarization (CALIOP). Since then, CALIPSO has been collecting an almost continuous 
record of high-resolution (333 m in the horizontal and as fine as 30 m in the vertical) profiles of 
aerosols and clouds attenuated backscatter at 532 and 1064 nm wavelengths, along with polarized 
backscatter in the visible channel, day and night, covering 82°N and 82°S (Winker et al., 2009; 2010). 
The polarization measurement allows for the separation of non-spherical dust particles from other 
types of aerosol. CALIOP algorithms use all the information: backscatter, the color ratio (ratio of 532 
to 1064 nm) and depolarization signal to first distinguish between clouds and aerosols in the column 
(Liu et al., 2009), and then to identify different aerosol types: namely smoke, polluted continental, 
polluted dust, dust, clean continental, and clean marine, with respective extinction-to-backscatter ratio 
or lidar ratio (S) of 70, 70, 65, 40, 35, and 20 sr at 532 nm (Omar et al., 2009). Once aerosol type is 
detennined, aerosol extinction profiles are retrieved by applying the extinction-to-backscatter ratio for 
the particular aerosol type identified. Alternative algorithms have been explored that attempt to 
provide a direct route to retrieving aerosol extinction without first assuming an extinction-to- 
backscatter ratio. Lor example, using the reflected laser signal from cloud tops as a light source, the 
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extinction of absorbing aerosol overlaying a cloud can be retrieved directly (Hu et al., 2007; Chand et 
al, 2008). 

The major limitation of CALIOP measurement is its near-zero swath. However, successive profiles 
observed over a period day and night have played important role in characterizing aerosol transport 
events. Aggregation of observations over seasonal and yearly time scales yields a climatology of 
aerosol vertical distribution, which has been used to determine aerosol transport height and provide 
vertically resolved mass fluxes. 

B. A movie: Aerosol intercontinental transport in 2001 as captured by MODIS and 
GOCART 

The movie illustrates the day-to-day evolution of aerosol optical depth and fractional contributions of 
fine and coarse particles for 2001. Top is derived from MODIS aerosol retrievals. Each frame 
represents a 5-day Gaussian weighted average centered on the indicated day. Bottom is derived from 
GOCART model simulations. Green indicates coarse mode dominated aerosol (dust), while red 
indicates fine mode dominated aerosol (smoke and pollution). Intensity of green and red increases as 
AOD increases. Individual images and movie were created by Reto Stockli under the umbrella of 
NASA's Earth Observatory and with input from Yoram Kaufman. 
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